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ABSTRACT 


-V 

This  thesis  is  a  study  of  performance  measures  for 
acoustooptic  spectrum  analysis,  and  optical  excision  of 
narrowband  interference.  A  Gaussian  laser  beam  model  is 
developed  and  programmed  for  propagation  through  thin  lenses. 
Effects  of  beam  truncation  by  Bragg  cell  apertures  are  ex¬ 
amined.  Bragg  cell  performance  is  analyzed  with  respect  to 
diffraction  efficiency,  information  capacity,  and  response 
agility.  A  performance  comparison  is  made  between  charge- 
coupled  device  and  photodiode  technologies  in  photodetector 
array  applications.  Overall  excisor  performance  is  discussed 
in  terms  of  interference  removal  effectiveness  and  process 
corruption  of  residual  information  bearing  signals. 
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I .  INTRODUCTION 


A.  PROBLEM  OVERVIEW 

The  increased  use  of  spread  spectrum  communication  and 
ranging  techniques  has  necessitated  rethinking  of  intercep¬ 
tion  scheme  implementations.  Two  common  requirements  of  many 
receiver  systems  operating  in  this  wideband  environment  are: 
adequate  spectrum  analysis,  and  effective  excision  of  narrow- 
band  high  power  interference.  Acoustooptic  spectrum  analysis 
and  optical  interference  excision  are  two  demonstrable  methods 
of  fulfilling  these  respective  requirements. 

An  acoustooptic  spectrum  analyzer  transforms  a  wideband 
time  domain  signal  into  a  spatially  detected  frequency  domain 
signal.  An  optical  excisor  at  the  detector  can  also  utilize 
this  space-time  transformation.  It  performs  real-time  notch 
filtering  of  narrowband  interference.  Specific  operational 
applications  of  the  two  devices  include  broadband  receiver, 
radar  and  surveillance  systems.  All  three  military  services 
rely  on  these  systems  and  need  good  interference  rejection. 

To  meet  this  need,  a  NAVELEX  350/DARPA  program  was  funded  for 
the  development  of  optical  excision  techniques. 

The  technical  objectives  of  this  effort  are  to: 

1.  Survey  potential  applications  and  analyze  technical 
requirements  imposed  by  each. 

2.  Analyze  alternative  optical  excision  configurations 
and  components . 
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3.  Develop/procure  major  system  components. 

4.  Fabricate  and  measure  alternative  broadband  excisor 
configurations  to  verify  analytical  results. 

5.  Fabricate  and  field  test  a  brassboard  model. 

6.  Design  a  preliminary  advanced  development  model. 

The  principal  civilian  contractor  for  the  project  is  PROBE 
Systems,  Inc.  PROBE,  assisted  by  the  development  efforts  at 
the  Naval  Research  Laboratory  (NRL)  and  Naval  Postgraduate 
School  (NPS) ,  was  tasked  with  attainment  of  these  objectives. 
NRL's  specific  program  assistance  was  the  development  of  a 
photodichroic  optical  clipper  to  be  evaluated  by  PROBE.  The 
research  summarized  in  this  thesis  is  part  of  the  NPS  con¬ 
tribution  to  the  development  program. 

B.  THESIS  OBJECTIVE 

The  objective  of  this  work  effort  was  to  conduct  a  per¬ 
formance  analysis  of  acoustooptic  (a-o)  spectrum  analyzer 
and  optical  excisor  components.  The  choice  of  this  goal, 
over  other  possible  goals  related  to  the  development  program, 
allows  maximal  use  of  the  Postgraduate  School's  substantial 
technical  library  and  computer  facilities.  The  approach  was 
also  felt  to  be  the  best  way  to  augment  PROBE'S  experimental 
study  and  minimize  duplicate  effort.  This  resulting  study 
is  intended  to  benefit  the  program  by  providing  a  comprehen¬ 
sive  set  of  "first  cut"  performance  guidelines. 

This  objective  was  divided  into  four  specific  aims.  Each 
effort  relates  to  a  particular  system  component.  A  summary 
of  these  goals  is  provided  in  Section  E.  However,  more 
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complete  explanations  of  a-o  spectrum  analysis  and  optical 
excision  are  required  beforehand. 

C.  ACOUSTOOPTIC  SPECTRUM  ANALYSIS 

An  acoustooptic  spectrum  analyzer  is  primarily  composed 
of  an  a-o  beam  deflector,  photodetector  array,  and  a  laser. 

Its  operating  sequence  begins  with  an  input  rf  signal  which 
excites  a  piezoelectric  transducer  mounted  on  the  acoustic 
medium  of  the  a-o  deflector  or  Bragg  cell.  The  transducer 
launches  an  acoustic  wave train  through  the  medium.  The  wave- 
train  spatially  modulates  the  refractive  index  and  acts  as  a 
moving  diffraction  grating  for  the  laser  light.  The  acousto¬ 
optic  interaction  produces  a  one-dimensional  angular  pattern 
of  diffracted  beams.  This  pattern  matches  the  spectral  con¬ 
tent  of  the  original  rf  signal.  Each  individual  beam  in¬ 
tensity  is  directly  proportional  to  the  corresponding  rf 
signal  frequency  component  intensity.  Beam  deflection  angles 
are  directly  proportional  to  the  component  frequencies.  A 
thin  convex  lens  is  used  to  focus  the  diffraction  pattern  on 
an  image  plane  linear  detector  array.  Lastly,  the  array 
photosites  video-detect  the  beam  pattern  and  transfer  the 
spectral  information  out  electronically  via  shift  registers. 

This  method  of  wideband  spectrum  analysis  offers  many 
features  which  are  difficult  to  achieve  by  other  means.  For 
example,  a  conventional  scanning  spectrum  analyzer  suffers 
from  several  shortcomings  when  used  as  a  wideband  surveillance 
receiver.  Signals,  which  utilize  direct  sequence  and  frequency 
hopping  spread  spectrum  modulation  techniques,  prove  difficult 
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to  adequately  monitor.  Pulsed  signals  are  also  troublesome 
because  their  probability  of  intercept  is  low.  An  a-o  spec¬ 
trum  analyzer  readily  accomplishes  these  surveillance  tasks. 
Since  the  instantaneous  bandwidth  is  very  wide,  probability 
of  intercept  problems  are  greatly  diminished.  The  demonstra¬ 
ted  output  frequency  resolution  is  excellent.  These  two 
features  permit  highly  effective  operation  in  densely  popu¬ 
lated  signal  environments.  Also,  the  simpler  physical  com¬ 
ponents  make  the  resulting  equipment  package  smaller  and  less 
expensive  than  comparable  conventional  systems,  especially  in 
the  case  of  integrated  optic  system  versions. 

Several  articles  which  discuss  currently  available  com¬ 
mercial  units  are  listed  in  the  bibliography.  Reference  1 
describes  an  ITEK  a-o  spectrum  analyzer  with  a  500  MHz  band¬ 
width  and  1  MHz  frequency  resolution.  Typical  applications 
for  this  processor  include  broadband  analysis  and  cueing 
receivers . 

Reference  2  discusses  GTE  Sylvania's  work  to  develop 
acoustooptic  wideband  receiver  and  direction  finder  systems, 
presenting  an  angle  of  arrival  detection  scheme  which  utilizes 
multichannel  Bragg  cells.  A  GTE  receiver  having  a  1  GHz  band¬ 
width  and  frequency  resolution  of  2.6  MHz  is  also  mentioned. 

A  snapshot  overview  of  acoustooptic  development  programs 
is  related  in  Ref.  3.  The  article  discusses  recent  advances 
made  by  the  above  two  companies  and  related  work  done  by 
others.  It  also  examines  exploratory  work  in  integrated  optic 
circuit  chips.  Hughes,  Rockwell  International,  and  Westing- 
house  are  the  major  companies  deeply  involved  in  this. 
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In  addition  to  the  above  articles,  the  bibliography  also 
notes  a  few  landmark  research  papers.  Reference  4  is  the 
first  article  in  this  category.  It  is  an  excellent  review 
of  acoustooptic  historical  aspects  which  covers  progress 
beginning  with  Brillouin's  1922  prediction  of  light  scattering 
by  thermally  generated  acoustic  waves.  In  this  work,  C.  G. 
Quate  and  his  associates  neatly  condensed  much  of  the  primarily 
academic  work  done  prior  to  the  1960 's.  They  also  conducted 
a  thorough  two-dimensional  analysis  of  wave  interactions. 

The  invention  of  the  laser  in  the  early  sixties  signifi¬ 
cantly  increased  the  practical  application  possibilities  of 
acoustooptic  devices.  Some  of  the  first  attempts  to  utilize 
the  Bragg  cell  are  discussed  in  Ref.  5.  The  authors,  A. 

Korpel  and  others,  explored  the  feasibility  of  acoustooptic 
television  displays.  They  presented,  in  addition  to  the 
experimental  findings,  an  analysis  of  one  of  the  first  a-o 
beam  deflectors  having  a  stepped  transducer.  The  special 
transducer  steered  the  acoustic  wavetrain  to  maximize  the 
proper  light-sound  interaction.  This  "beam  steering"  tech¬ 
nique  increased  the  effective  operating  bandwidth  of  the 
Bragg  C611. 

R.  Adler,  one  of  the  pioneers  who  worked  with  Korpel, 
gave  a  clear  summary  of  beam  deflector  principles  in  Ref.  6. 

He  also  discussed  a  wide  variety  of  acoustooptic  applications 
which  included  television  displays,  spectrum  analyzers,  and 
acoustic  imagers. 

Other  references,  which  relate  to  specific  analyzer 
aspects,  are  cited  in  following  sections.  However,  the  above 
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three  papers  best  provide  the  flavor  of  development  work  and 
concisely  present  the  significant  acoustooptic  principles. 

D.  OPTICAL  EXCISION 

The  optical  excision  process  utilizes  all  a-o  spectrum 
analyzer  components  with  the  addition  of  an  optical  or  elec¬ 
tronic  clipper.  Only  minor  modifications  are  made  in  the 
operating  sequence.  The  image  plane  diffraction  pattern  is 
generated  in  the  same  manner  as  previously  discussed.  The 
clipping  action  takes  place  in,  or  in  front  of,  the  image 
plane.  Phase  integrity  is  important  for  the  time  domain 
reconstruction  of  the  signal,  thus  heterodyne-detection 
replaces  video-detection  to  preserve  signal  phase. 

If  an  electronic  clipper  is  used,  then  the  array  photo¬ 
site  outputs  must  be  individually  serviced.  Outputs  which 
correspond  to  frequencies  in  the  notch  band  are  ignored.  The 
remaining  outputs  are  summed  to  produce  an  excised  version 
of  the  input  signal. 

An  optical  clipper  requires  a  single  large  area  photo¬ 
detector  to  recover  the  heterodyned  signal.  An  opaque  ob¬ 
struction  blocks  or  clips  the  diffracted  beams  produced  by 
frequencies  in  the  notch  band.  The  remaininq  beams  are 
mixed  with  a  reference  local  oscillator  beam  and  recovered 
by  the  detector.  The  net  output  is  again  an  excised  version 
of  the  input  signal.  A  diagram  of  the  optical  clipping  pro¬ 
cess  is  shown  in  Figure  1. 

Optical  excision  has  several  advantages  lacking  in 
conventional  approaches  to  narrowband  interference  rejection. 


Optical  Clipping  Process  Sequence 


The  process  offers  an  arbitrary  notch  bandwidth  with  nearly 
infinite  skirt  slope.  Optical  excision  lends  itself  well  to 
adaptive  notch  filtering,  which  is  important  in  signal  environ¬ 
ments  having  time-varying  interference  sources  of  arbitrary 
number  and  center  frequencies. 

Several  efforts  to  develop  suitable  clippers  are  under¬ 
way.  Photodichroic  crystals  are  one  set  of  possible  con¬ 
tenders  for  application  in  optical  clippers.  These  crystals, 
when  used  with  crossed  polarizers,  act  as  optical  attenuators. 
The  amount  of  attenuation  is  dependent  upon  the  incident  light 
intensity.  This  property  gives  them  an  inherent  spatial  fil¬ 
ter  capability.  Strong  interference  signals  are  sharply 
attenuated  or  clipped,  and  weaker  signals  of  interest  pass 
through.  As  mentioned  in  the  overview,  NRL  is  deeply  invol¬ 
ved  in  this  development  area.  Reference  7  provides  a  summary 
of  experimental  work  on  photodichroic  spatial  filters  per¬ 
formed  by  W.  C.  Collins  at  the  Naval  Research  Laboratory. 

PROM's  or  Pockels  Readout  Optical  Modulators  also  show 
promise  in  this  area.  A  PROM  optical  clipper  uses  two  lasers 
of  different  wavelengths.  The  first  or  "write"  laser  has  a 
shorter  wavelength  and  creates  photoconductive  spots  on  the 
PROM.  These  areas  correspond  to  the  interference  frequencies 
to  be  excised.  The  PROM,  positioned  between  crossed  polar¬ 
izers,  is  "read"  by  a  laser  of  longer  wavelength.  Diffracted 
beams  produced  by  the  interference  signal  are  blocked  by  the 
combined  action  of  the  photoconductive  regions  and  crossed 
polarizers.  The  remaining  diffracted  beams  are  transmitted 
through  the  PROM  and  recovered  by  the  heterodyned  detector. 
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Reference  8  describes  a  PROM  developed  by  ITEK  which  was  used 
as  part  of  an  acoustooptic  spectrum  analyzer. 

The  LCLV  or  Liquid  Crystal  Light  Valve  is  another  possible 
optical  clipper.  Hughes  Research  Laboratory  has  been  the 
principal  developer.  The  device  consists  of  a  liquid  crystal 
photoconductor  sandwich  which  uses  scattering  centers  to  dis¬ 
rupt  the  coherence  of  incident  laser  light.  The  scattering 
centers  are  created  by  the  interaction  of  a  cross-sandwich 
electric  field  and  the  incident  light  intensity  pattern.  The 
incident  light  lowers  the  resistivity  of  the  photoconductor 
sandwich  portion  in  proportion  to  the  intensity.  This  action 
increases  the  electric  field  in  the  liquid  crystal  layer  and 
causes  molecular  alignment  rotations.  The  areas  of  molecular 
rotation  are  the  scattering  centers.  In  order  to  be  used  as 
an  optical  clipper,  the  light  valve  is  positioned  between 
crossed  polarizers.  Attenuation  is  dependent  upon  molecular 
rotation  and  thereby  the  incident  light  intensity.  Reference 
9  discusses  progress  made  by  Hughes  on  the  LCLV.  The  article 
also  compares  various  spatial  filters'  performance  including 
the  three  mentioned  here. 

E.  SPECIFIC  GOALS 

The  specific  goals  of  this  work  were  narrowed  down  to 
four  principal  areas.  First,  laser  beam  profile  effects  on 
signal  frequency  component  resolution  required  study.  Second, 
Bragg  cell  performance  factors  needed  determination  for  optimal 
system  design.  Third,  photodetector  capabilities  and  effects 
required  scrutiny  to  compare  the  advantages  of  competing 
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detector  technologies.  And  fourth,  optical  excision  needed 
modeling  to  determine  process  effectiveness.  The  following 
paragraphs  expand  on  these  points . 

Laser  light  beam  propagation  was  the  first  major  aspect 
studied.  This  subject  is  important  because  the  beam's  cross- 
sectional  behavior  is  a  function  of  the  propagation  circum¬ 
stances.  The  beam  profile  significantly  affects  both  a-o 
spectrum  analyzer  and  optical  excisor  performance. 

For  example,  a  wide  beam  is  desired  at  the  deflector 
aperture.  This  condition  maximizes  the  acoustooptic  inter¬ 
action.  Conversely,  a  small  light  beam  cross  section  is 
necessary  at  the  image  plane  to  discriminate  the  input 
signal's  frequency  components. 

Ray  optics  can  not  be  used  to  determine  beam  profile 
behavior  through  the  various  system  elements,  as  it  does  not 
adequately  describe  beam  propagation  characteristics.  For¬ 
tunately,  the  complex  variable  equations,  which  do  govern 
beam  transmission  through  space  and  thin  lenses,  are  avail¬ 
able;  however,  they  are  cumbersome  to  use.  Simple  methods 
of  dealing  with  them  were  needed  to  increase  their  useful¬ 
ness  as  performance  gauges.  In  particular ,  a  hand-held  cal¬ 
culator  program  was  desirable  to  determine  optimal  lens 
placement  and  beam  propagation  indices. 

Also,  propagation  through  the  Bragg  cell's  finite 
dimension  aperture  affects  cross-sectional  behavior.  The 
abrupt  windowing  action  truncates  the  beam  profile  and 
causes  sidelobes  to  appear.  These  sidelobes,  if  sufficiently 
large,  contribute  erroneous  frequency  components  to  the  image 
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plane.  Investigation  of  this  phenomenon  was  required  to 
determine  the  degree  of  performance  degradation. 

The  second  thesis  goal  was  to  analyze  factors  affecting 
Bragg  Cell  performance.  To  pursue  this  goal,  it  was  neces¬ 
sary  to  define  good  deflector  performance.  Uchida  /— 10 _7 
did  this  concisely  in  terms  of  three  criteria:  high 
diffraction  efficiency,  large  information  capacity,  and  quick 
response . 

Diffraction  efficiency  is  the  ratio  of  a  diffracted  light 
beam's  intensity  to  that  of  the  incident  light  beam.  The 
input  acoustic  power  is  one  of  the  components  which  determine 
diffraction  efficiency.  A  higher  input  acoustic  power  will 
give  a  greater  diffraction  efficiency,  but  only  at  a  price. 
Increased  acoustic  power  heats  the  interaction  medium  and 
creates  troublesome  temperature  gradients.  Resolution 
suffers  and  information  capacity  is  lowered.  Damping  of 
internal  reflections  is  made  more  difficult,  causing  an  un¬ 
desirable  persistence  of  signals  which  are  no  longer  present. 
Deflectors  which  exhibit  high  diffraction  efficiency  with 
low  acoustic  power  are  very  advantageous. 

Information  capacity  can  be  measured  simply  as  the  number 
of  resolvable  diffracted  beams.  It  can  also  be  influenced 
by  geometric  factors.  Stretching  the  interaction  medium's 
acoustic  path  length  will  increase  the  number  of  resolvable 
frequencies.  However,  this  approac.i  will  also  slow  respon¬ 
siveness  to  changing  signal  environments. 

Simultaneous  achievement  of  all  three  criteria,  high 
diffraction  efficiency,  large  information  capacity,  and  quick 
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response  is  not  possible  yet,  but  performance  may  be  adequately 
tailored  to  a  particular  application  by  decreasing  one  un¬ 
necessary  feature  in  order  to  increase  a  crucial  one. 

Photodetector  performance  was  another  important  subject 
for  analysis.  Detector  performance  criteria  include:  dynamic 
range,  resolution,  and  information  bandwidth.  For  spectrum 
analysis  applications,  there  are  two  contending  device  types. 
Either  a  Charge-Coupled  Device  (CCD)  detector  array  or  a 
Charge-Coupled  PhotoDiode  (CCPD)  detector  array  may  be  used. 
Both  arrays  are  composed  of  a  linear  grouping  of  photosites. 

The  photosites  detect  and  integrate  incident  light  intensities. 
Their  response  sensitivity  is  matched  to  the  laser  wavelength. 
Dynamic  range  is  limited  by  photosite  saturation  illumination 
and  dark  current.  Resolution  is  limited  by  the  number  and 
linear  packing  density  of  the  photosites.  Information  band¬ 
width  is  limited  by  electronic  sampling  considerations. 

For  optical  excision  applications  which  utilize  an  opti¬ 
cal  clipper,  a  single  large  area  detector  is  required.  Po¬ 
tential  device  types  include  photodiodes  and  photomultiplier 
tubes.  A  photodiode  may  be  adequate  for  a  particular  system 
if  the  information  bandwidth  of  the  signal  is  not  large. 
Photomultiplier  tubes  are  available  for  wideband  requirements. 
Dynamic  range  limitations  are  similar  to  those  for  detector 
arrays.  Resolution  is  determined  by  the  optical  clipper  and 
the  diffracted  beam  spot  size. 

A  number  of  references  are  available  which  highlight 
important  aspects  of  the  different  detectors.  CCD  principles 
of  operation  are  presented  in  Ref.  11.  Reference  12  analyzes 
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photodiode  array  noise  considerations  in  depth.  Reference 
13  includes  mathematical  examinations  of  video  and  heterodyne- 
detection  for  photodiodes  and  photomultiplier  tubes. 

Electronic  and  optical  clipping  methods  comprised  the 
final  areas  of  analysis.  Models  of  both  techniques  were 
developed.  The  first  model  addressed  optical  clipping. 

The  blocking  action  of  an  optical  clipper  may  be  modeled 
to  a  first  order  approximation  as  a  narrow  obstacle  in  a 
coherent  plane  wave  light  field.  The  corresponding  image 
plane  diffraction  pattern  can  be  described  in  terms  of 
Fresnel  integrals.  These  integrals  may  be  approximated  by 
numerical  methods.  The  resulting  computer  plots  help  one  to 
visualize  diffraction  effects  produced  by  an  optical  clipper. 

The  next  model  is  primarily  associated  with  electronic 
clipping.  It  is  intended  to  provide  a  measure  of  signal 
degradation  which  is  induced  by  the  excision  process.  This 
degradation  which  is  caused  by  the  removal  of  the  desired 
signal's  frequency  components  in  the  excisor's  notch  band. 

For  Direct  Sequence  Spread  Spectrum  (DSSS)  signals,  this 
results  in  attenuation  of  the  correlation  peak.  By  compari¬ 
son,  degradation  effects  for  frequency  hopped  signals  are  not 
considered  here  as  these  signals  do  not  derive  any  benefits 
from  the  excision  process. 

F.  PERSPECTIVE 

The  fundamental  purpose  of  this  introduction  was  to  ex¬ 
plain  the  motivations  for  the  work  effort.  Primary  justifi¬ 
cation  for  the  undertaking  was  the  same  as  for  the  sponsoring 
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development  program — a  real  obstacle  to  effective  utilization 
of  spread  spectrum  techniques  must  be  confronted.  The  pro¬ 
blem  of  narrowband  interference  cannot  reasonably  be  expected 
to  diminish  with  only  passage  of  time,  more  probably  it  will 
increase . 

Acoustooptics  promises  a  solution  to  this  interference 
problem.  The  necessary  technologies  are  maturing,  and  several 
available  products  ably  handle  missions  similar  to  the  in¬ 
tended  ones  discussed  here.  Therefore,  this  promised  solu¬ 
tion  is  very  likely  to  be  fulfilled  in  the  near  future. 
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II.  PERFORMANCE  GUIDELINE  ANALYSIS 


A.  PARAMETERS  OF  NEGLIGIBLE  EFFECT 

Several  secondary  effects  are  present  in  both  the  a-o 
spectrum  analyzer  and  optical  excisor  which  were  not  con¬ 
sidered  as  being  germane  to  the  thesis.  Consequently,  they 
were  not  examined  in  great  detail.  This  section  is  a  brief 
compilation  of  those  effects. 

One  performance  affecting  parameter  outside  the  analysis 
realm  of  this  thesis  is  input  noise  contributed  by  the  re¬ 
ceiver  system.  No  operation  in  either  spectrum  analysis  or 
optical  excision  processors  decreases  white  Gaussian  noise. 
Above  50  MHz,  atmospheric  noise  diminishes,  and  receiver  noise 
figure  becomes  an  important  consideration.  The  intended  appli¬ 
cations  for  both  systems  are  assumed  to  be  above  50  MHz.  How¬ 
ever,  the  techniques  dealing  with  receiver  noise  figure  are 
well  documented  elsewhere  and  therefore  are  not  addressed  in 
this  report.  In  all  further  analysis,  the  system  noise  in¬ 
duced  prior  to  Bragg  cell  injection  is  assumed  negligible. 

External  vibration  is  another  potential  problem  in  both 
a-o  spectrum  analyzers  and  optical  excisors.  Component 
mounting  vibrations  can  introduce  undesirable  signal  modula¬ 
tions.  Air  currents  may  also  disturb  the  signal.  Coherent 
detection  schemes  are  particularly  sensitive  to  both  occur¬ 
rences  .  Susceptibility  to  small  ambient  air  currents  was 
documented  by  M.  King  and  his  associates  in  Ref.  14.  They 
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observed  random  phase  modulation  of  a  coherently  detected 
signal  which  was  caused  by  minor  air  movement.  The  experi¬ 
mental  work  of  this  thesis  was  done  under  controlled  condi¬ 
tions  with  minimal  mechanical  disturbance.  Also,  video¬ 
detection,  which  reduced  vibration  concerns,  was  employed. 

In  general,  vibration  considerations  cannot  be  effectively 
encompassed  within  the  framework  of  this  report.  The  subject 
must  be  addressed  with  respect  to  each  particular  application. 

Laser  frequency  stabilization  requirements  of  the  optical 
excisor  coherent  detection  scheme  were  also  briefly  investi¬ 
gated.  Stabilization  was  concluded  to  be  unnecessary  for  this 
application  as  the  reference  local  oscillator  is  derived  from 
the  same  laser  beam  as  the  information  signal.  Frequency 
variations  in  the  local  oscillator  are  matched  by  those  in 
the  signal  beam.  Only  if  there  were  a  large  path  length  dif¬ 
ference  between  the  signal  beam  and  local  oscillator  would 
there  be  any  chance  of  beat  frequencies  appearing  in  the 
information  band.  This  condition  does  not  exist  in  the  opti¬ 
cal  excisors  previously  described. 

Tyndall  scattering  is  a  consideration  in  systems  which 
employ  liquid  medium  Bragg  cells.  It  is  not  a  concern  with 
high  purity  crystal  medium  Bragg  cells.  Tyndall  scattering 
is  due  to  foreign  suspended  matter  in  the  interaction  medium. 
Scattering  of  this  type  contributes  to  background  noise.  It 
can  be  controlled  by  careful  filtering  of  the  liquid  medium. 

Liquid  medium  Bragg  cells  are  restricted  to  operation 
under  100  MHz  due  to  attenuation  considerations.  Most  new 
Bragg  cells  employ  low  acoustic  velocity  crystals.  The  a-o 
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deflector,  used  in  experiments  related  to  this  thesis,  housed 
a  crystalline  interaction  medium,  therefore,  discussion  of 
this  topic  is  limited  here.  For  these  interested,  further 
remarks  on  Tyndall  scattering  are  made  in  Ref.  15. 

Temperature  gradients  in  the  Bragg  cell  medium  have  al¬ 
ready  been  described  in  the  introduction  as  harmful  to  infor¬ 
mation  capacity.  Reference  10  documents  experimental  obser¬ 
vations  of  thermal  gradient  effects.  This  condition  can  be 
controlled  in  essentially  three  ways.  First,  keeping  the 
input  acoustic  power  low  minimizes  input  heat.  Second,  em¬ 
ploying  materials  with  high  thermal  conductivity  helps  to 
maintain  uniform  heat  content  throughout  the  medium.  Third, 
using  external  cooling  aids  in  dissipating  heat  introduced 
by  the  transducer. 

Nonuniformity  of  photosite  responsivity  for  photodetector 
arrays  is  the  last  topic  considered  in  tnis  section.  Typical 
variation  for  both  CCD  and  CCPD  arrays  is  ten  percent  above 
or  below  the  mean.  Normally,  the  intended  applications  do 
not  require  more  exacting  tolerances.  In  the  event  that  such 
precision  is  required,  digital  compensation  techniques  may  be 
employed  to  offset  individual  photosite  deviations.  Reference 
12  contains  more  discussion  on  this  subject. 

B.  LASER  BEAM  PROPAGATION 

The  first  goal  of  the  propagation  study  was  to  develop 
simple  methods  of  dealing  with  the  equations  which  govern 
transmission  through  space  and  thin  lenses,  and  enable  deter¬ 
mination  of  beam  spot  size  and  phase  front  radius  anywhere 
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within  the  optic  train  path.  In  order  to  explain  the  results 
of  the  effort  more  concisely,  a  brief  review  of  laser  propa¬ 
gation  fundamentals  is  presented.  This  review  is  based  on  a 
treatment  provided  in  Ref.  16. 

The  type  of  laser  beam  source  discussed  throughout  this 
paper  is  a  standard  Helium-Neon  unit  which  radiates  at  a 
.6328  micron  wavelength.  The  beam  has  a  Gaussian  intensity 
profile  and  a  spherical  phase  front.  Beam  spot  size  is  de¬ 
fined  as  the  perpendicular  distance  from  the  beam  center- 

2 

line  to  a  point  where  the  intensity  has  dropped  by  1/e  factor. 
Beam  radius  is  the  phase  front  curvature  radius.  These  two 
quantities  uniquely  determine  the  characteristics  of  free 
space  Gaussian  beam  propagation.  The  plane  perpendicular  to 
the  propagation  path,  where  the  spot  size  is  a  minimum  and 
the  radius  is  infinite,  is  referred  to  as  the  "waist"  and  is 
shown  in  Figure  2. 

Two  equations  determine  beam  spot  size  and  radius  as  a 
function  of  directed  propagation  distance.  Plots  of  beam 
radius  and  spot  size,  for  propagation  through  free  space,  are 
depicted  in  Figs.  3  and  4.  The  equations  are: 

2  2 

Beam  Radius  =  R(z)  =  z  +  {Trwo/X}  /z  (1) 

Beam  Spot  Size  =  w(z)  =  w  {1  +  {Xz/{irw^}  }/  **  (2) 

c  o  o 


where : 

wq  =  Minimum  Beam  Spot  Size 

X  =  Laser  Wavelength 

z  =  Propagation  Distance  from  Waist 


Figure  2.  Gaussian  Beam  Characteristics 
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Figure  3. 
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Figure  4 .  Beam  Spot  Size  as  a  Function 
of  Propagation  Distance 
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In  order  to  predict  beam  behavior  through  a  thin  lens , 
the  quantity  q  must  be  introduced.  This  value  is  known  as 
the  'complex  beam  radius'.  It  is  defined  by  the  following 
relation. 

2 

Complex  Beam  Radius  =  l/q(z)  =  1/R(z)  -  jX/irw  (z)  (3) 

At  the  waist,  the  beam  radius  R(0)  goes  to  infinity,  leaving: 

q(0)  =  j7TW^/\  =  qQ  (4) 

Patient  manipulation  of  Equations  (1)  through  (4)  will  yield 
Equation  (5),  which  completely  describes  Gaussian  beam  trans¬ 
mission  in  free  space. 

q(z)  =  qQ  +  z  (5) 

Spherical  waves  which  pass  through  thin  lenses  obey  the 
lens  law.  For  Gaussian  beams,  the  complex  beam  radius  under¬ 
goes  the  following  transformation  shown  in  Figure  5. 

l/q2  =  1/q-L  -  1/f  (6) 

where : 

q^  =  incident  beam's  complex  radius 

q^  =  exit  beam's  complex  radius 

f  =  thin  lens  focal  length 

The  above  equations  provided  the  basis  for  the  TI-58  or 
TI-59  hand-calculator  program  to  deal  with  Gaussian  beam 
propagation  through  free  space  and  thin  lenses.  This  program 
is  listed  and  explained  in  Appendix  B.  The  first  practical 
application  of  the  program  was  to  compute  the  beam  spot  size 
incident  to  the  image  plane  photodetector  array  of  an 
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experimental  acoustooptic  spectrum  analyzer  constructed  as 
part  of  this  work  effort. 

The  calculation's  accuracy  was  dependent  upon  one  pri¬ 
mary  assumption.  The  beam  profile  could  not  be  significantly 
truncated  by  the  Bragg  cell  aperture  or  any  of  the  thin  lenses 
in  the  beam  path.  This  assumption  was  adequately  justified 
for  the  lenses  used.  A  lens  diameter  of  three  times  the  spot 
size  is  sufficient  to  contain  99  percent  of  the  beam  power. 

All  lenses  in  the  experimental  set-up  were  greater  than  five 
spot  sizes  wide. 

However,  this  criterion  was  not  upheld  as  the  beam  passed 
through  the  3ragg  cell  aperture.  The  beam  spot  size  was 
expanded  to  completely  illuminate  the  aperture  and  thereby 
enhance  information  capacity.  Consequently,  the  beam  pro¬ 
file  was  substantially  truncated.  In  the  limit,  the  incident 
beam  cross-sectional  intensity  would  have  become  uniform  in 
distribution.  The  resulting  exit  beam,  were  it  produced  by 
the  interaction  of  a  plane  sound  wave  and  a  rectangular  light 
aperture  of  finite  dimensions,  would  have  had  a  horizontal 
intensity  profile  given  below  by  Equation  (7) .  This  formula 
is  based  on  the  Fraunhofer  approximation  and  is  derived  in 
Ref.  17. 


I(r, 9) 


I  (r); 
o 


;sinf  XD9/X }  i  ^ 
1  {  ttD6 /X  ; 


(7) 


where : 

I(r,9)  =  Intensity  as  a  function  of  angle  off  the  Bragg 
cell  aperture  boresight  and  path  length  to  the 
image  plane 


32 


IQ(r)  =  Boresight  light  beam  intensity 
D  =  Aperture  horizontal  dimension 

X  =  Laser  wavelength 

9  =  Angle  off  boresight 

r  =  Distance  to  transform  lens 

2 

The  1/e  beamwidth  was  determined  from  Equation (7)  by  an 
iterative  algorithm.  The  solution  is  given  by  Equation  (8). 

2 

Sine  Function  1/e  Beamwidth  =  SBW  =  1.4r\/D  (8) 

In  order  to  determine  the  seriousness  of  violating  this 
criterion,  beamwidth  calculations  were  made  for  both  extremes. 
The  first  assumed  an  incident  Gaussian  beam  with  no  truncation. 
The  second  assumed  an  incident  plane  wave.  Results  of  the 
calculations  were  compared  with  an  experimentally  determined 
value . 

The  equipment  used  to  make  confirming  observations  of 
image  plane  beam  spot  size  is  block  diagrammed  in  Figure  6. 
Component  devices  are  described  in  the  following  paragraphs, 
and  dimensional  quantities  needed  by  the  propagation  program 
are  provided  in  Table  I. 

The  laser  utilized  for  the  measurements  was  a  Spectra- 
Physics  115.  It  produced  a  .5  milliwatt  continuous  wave 
output  which  had  a  TEMqo  spatial  mode  and  random  polarization. 
Power  line  ripple  and  random  noise  contributions  were  limited 
to  one  quarter  percent  of  the  output  power.  The  longitudinal 
mode  spacing  was  550  MHz. 
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Figure  6.  Experimental  Acoustooptic  Spectrum  Analyzer 
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TABLE  I.  PROPAGATION  PROGRAM 

LASER 

EXIT  BEAM  SPOT  SIZE 
EXIT  BEAM  RADIUS 
WAVELENGTH 

SPATIAL  FILTER 
DISTANCE  FROM  LASER 
LENS  FOCAL  LENGTH 

COLLIMATING  LENS 
FOCAL  LENGTH 

TPANSFORM  LENS 
FOCAL  LENGTH 

DISTANCE  FROM  COLLIMATING 
LENS 


INPUTS 

.44  mm 
INFINITY 
.6328  urn 

.  1  m 
14.8  mm 

.254  m 

.  304  m 

.  3  m 
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A  spatial  filter  followed  the  laser  in  the  operating 
train  sequence.  It  removed  intensity  fluctuations  which 
were  due  to  light  scattering  from  the  beam  profile.  Dust 
on  the  laser  windows  and  output  mirrors  is  the  major  cause 
of  this  disturbance,  although  lens  imperfection  may  be  a  con¬ 
tributing  factor  as  well. 

The  spatial  filter  performed  its  cleaning  task  by  focus¬ 
ing  the  beam  through  a  small  pinhole.  This  focusing  action 
caused  the  scattered  light  to  spread  out  from  the  central 
beam.  Light  remaining  in  the  central  beam  had  a  smoothed 
Gaussian  profile.  The  pinhole  passed  the  major  portion  of 
the  central  beam  and  removed  the  scattered  light.  Slight 
truncation  of  the  beam  by  the  pinhole  produced  an  Airy  Disk 
diffraction  pattern.  This  phenomenon  is  discussed  by  Born 
and  Wolf  in  Ref.  18. 

The  outer  rings  of  the  Airy  pattern  were  removed  by  pass¬ 
ing  the  central  beam  through  a  second  aperture  or  iris.  The 
net  far  field  beam  profile  had  a  nearly  ideal  Gaussian  shape. 
This  beam  was  collimated  by  a  thin  lens.  Doing  so  ensured 
complete  illumination  of  the  succeeding  Bragg  cell  aperture. 

The  particular  Bragg  cell  employed  in  this  set  up  was  an 
Isomet  (1205-1-2) .  It  used  Lead  Molybdate  (PbMo04)  as  the 
interaction  medium.  The  bandpass  center  frequency  was  80 
MHz.  This  translated  to  a  deflection  angle  of  7  milliradians . 
Nominal  rf  input  power  was  one  and  a  half  Watts. 

An  80  MHz,  fixed  amplitude,  continuous  sine  wave  was  chosen 
as  the  input  drive  signal.  A  pure  sinusoidal  tone  was  used  so 
that  no  part  of  the  measured  beamwidth  was  due  to  the  input 
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signal.  The  deflected  beam  was  assumed  to  be  a  function  only 
of  the  interaction  between  the  Gaussian  light  beam  and  the 
Bragg  cell  aperture's  size  and  shape.  Discussion  of  profile 
affecting  light-sound  interplay  is  delayed  until  the  next 
section . 

When  the  Bragg  cell  was  excited  by  the  input  signal,  two 
light  beams  exited  the  device  orifice.  The  first  beam  or 
"zero  order  diffraction  beam"  was  not  deflected.  The  only 
information  carried  by  it  was  in  the  form  of  amplitude  modu¬ 
lation.  A  stronger  input  signal  would  have  meant  a  weaker 
undeflected  beam.  This  beam  was  not  useful  and  therefore 
terminated  by  a  opaque  zero  order  stop.  The  second  beam  was 
deflected,  and  represented  first  order  diffraction.  It  was 
passed  through  the  transform  lens  and  brought  to  a  focus  on 
the  image  plane  linear  detector  array  for  width  measurement. 

The  transform  lens  played  either  one  of  two  different 
roles  which  depended  upon  whether  the  diffracted  beam  had  a 
Gaussian  or  Sine  function  intensity  profile.  If  the  beam 
were  truly  Gaussian  and  still  well  collimated  as  it  exited 
the  Bragg  cell,  then  the  image  plane  spot  size  would  have 
been  insensitive  to  changes  in  the  distance  between  Bragg 
cell  and  the  transform  lens.  The  important  distance  which 
needed  precise  calculation  would  have  been  that  from  the  lens 
to  the  image  plane.  In  the  case  of  a  well  collimated  Gaussian 
beam,  this  span  is  nearly  equal  to  the  lens  focal  length. 

Conversely,  had  the  exit  beam  exhibited  a  Sine  function 
intensity  profile,  then  the  distance  from  the  lens  to  the 
image  plane  would  have  been  arbitrary.  Such  a  beam  would 
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have  appeared  to  emanate  from  a  point  source  at  the  Bragg 
cell  aperture.  A  thin  lens  positioned  a  focal  length  away 
would  have  collimated  the  beam.  The  resulting  spot  size, 
as  shown  in  Equation  (7) ,  would  have  been  a  function  only 
of  the  distance  to  the  transform  lens,  and  not  of  the  trans¬ 
form  lens-image  plane  separation. 

The  observed  image  plane  spot  size  was  sensitive  to 
changes  in  either  separation  distance.  This  result  implied 
the  exit  beam  had  retained  characteristics  of  both  profiles. 
The  optimal  transform  lens  position  was  located  approximately 
one  focal  length  away  from  the  Bragg  cell  and  a  focal  length 
from  the  image  plane.  This  "best  arrangement"  is  confirmed 
by  Whitman  and  his  associates  in  Ref.  19. 

The  image  plane  sensor  used  to  detect  beamwidth  was  a 
(CCD121H)  linear  detector  array  developed  by  Fairchild.  It 
consists  of  1728  sensor  elements  and  is  a  charge  coupled, 
buried  channel  device.  The  output  is  delivered  serially 
and,  for  the  experiment,  went  to  a  Tektronics  422  oscillo¬ 
scope  display.  A  typical  output  display,  which  resulted  from 
an  80  MHz  monotone  input  signal,  is  shown  in  Figure  7. 

The  physical  apportionment  of  the  photo  elements  is 

2 

shown  m  Figure  8 .  The  1/e  beamwidth  was  measured  by  ex¬ 
panding  the  oscilloscope  time  display  and  counting  photosite 
output  peaks  within  the  beamwidth.  This  number  was  converted 
into  a  spatial  dimension  by  using  the  known  dimensions  in 
Figure  8. 
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Figure  7.  Photodetector  Array  Output 
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The  Sine  profile  was  determined  using  Equation  (8)  and 
the  following  values: 

r  =  focal  length  of  transform  lens  =  .304  m 
X  =  Helium-Neon  Laser  wavelength  =  .6328  um 

D  =  effective  aperture  width  =  .7  cm 

The  Gaussian  beamwidth  was  calculated  by  using  the  propaga¬ 
tion  program  and  Table  I.  Intermediate  results  are  listed 
in  Table  II.  Incidently,  this  table  may  be  used  as  a 
standard  solution  set  to  ensure  all  proper  entries  are  keyed 
into  the  calculator's  program  memory. 

Table  III  compares  the  three  independently  determined 
values.  These  results  indicate  the  actual  beam  profile  was 
very  similar  to  the  Sine  function.  However,  the  highest 
sidelobe  was  measured  17  dB  down  from  the  central  peak.  If 
the  profile  were  an  ideal  Sine  function,  then  the  highest 
sidelobe  would  have  been  approximately  14  dB  down. 

D.  L.  Hecht  presented  a  concise  sidelobe  analysis  of 
truncated  Gaussian  beams  in  Ref.  20.  In  one  graph  he  related 
the  maximum  sidelobe  level  to  beam  truncation  ratio-  D/2w. 
This  graph  is  shown  in  Figure  9.  Again,  this  does  not  take 
into  account  effects  of  the  sound  column  of  the  light  beam 
profile,  it  is  a  function  only  of  truncation  ratio. 

Figure  9  provided  a  rough  means  of  indirectly  determining 
truncation  ratio.  For  a  -17  dB  sidelobe,  the  corresponding 
truncation  ratio  was  in  the  range  of  .55  to  .6.  The  actual 
truncation  ratio  value  was  calculated  as  .47.  It  was  deter¬ 
mined  by  using  the  spot  size  given  in  Table  II  for  the  beam 
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TABLE  II.  BEAM  PROFILE  INTERMEDIATE  RESULTS 

At  Laser  Aperture 

Beam  Spot  Size  =  w  =  .44  mm 

Beam  Radius  =  R  =  °° 

After  propagating  .1  m 

w  =  .4424  mm 

R  =  9.377  m 

At  Spatial  Filter  Lens  Exit 

■ 

w  =  .4424  mm 
R  =  -14 .82  mm 

After  propagating  14.820044  mm 

w  =  6.749  urn 

R  =  oo 

r 

\  ' 

After  propagating  .254  m 

4  w  =  7.581  mm 

R  =  .2540 

At  Collimating  Lens  Exit 

w  =  7.581  mm 
R  =  -3.205  x  IQ5  m 

After  propagating  .3  m 

w  -  7.581  mm 
R  =  1.770  x  106  m 

At  Transform  Lens  Exit 

w  =  7.581  mm 
R  =  -.3040  m 

.  After  propagating  3.039997  m 

w  =  8 . 077  urn 

R  =  oo 
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TABLE  III.  COMPARISON  OF  BEAMWIDTH  VALUES 


Sine  Image  Plane  Beamwidth  38  um 
Measured  Image  Plane  Beamwidth  34  um 
Gaussian  Image  Plane  Beamwidth  16  um 
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Figure  9.  Maximum  Sidelobe  Level  vs  Trancation  Ratio 


entering  the  Bragg  cell  (.75  cm)  and  the  deflector  aperture 
width  (.7  cm).  The  graph  determined  value  differed  by  17 
percent. 

This  disagreement  is  not  very  surprising  in  view  of  pos¬ 
sible  sound  column  effects  and  other  error  sources.  Channel 
stops  were  believed  to  have  played  a  role  in  limiting  pre¬ 
cision  measurements.  The  intervening  metal  strip  between 
two  adjacent  photosites  consumes  24  percent  of  the  available 
surface  area.  Its  effect  can  be  minimized  but  not  eliminated 
by  "rocking"  the  diffracted  beam  back  and  forth. 

Lack  of  a  precision  focus  could  have  also  distorted  ob¬ 
servations.  If  present,  minor  lens  and  Bragg  cell  imperfec¬ 
tions  would  have  slightly  affected  ultimate  focusability  of 
the  beam,  however,  operator  error  was  not  particularly  sus¬ 
pect  here.  Even  though  focusing  was  a  tedious  exercise, 
minor  lens  position  changes  produced  very  noticeable  amplitude 
deviations  in  an  expanded  scale  oscilloscope  display. 

Also,  the  detector  array  was  not  a  perfectly  one-dimen¬ 
sional  sensor.  Each  pixel  averaged  the  incident  light  in¬ 
tensity  over  a  rectangular  area.  The  vertical  photosite 
dimension  was  twice  the  magnitude  of  the  horizontal.  It  was 
also  equal  to  the  main  lobe  beamwidth  calculated  for  an  ideal 
Gaussian  beam.  It  is  not  fair  to  assume  that  the  ratio  of 
averaged  intensities  would  have  necessarily  been  identical  to 
the  ratio  of  the  peak  values.  A  closer  ratio  agreement  would 
have  been  realized  with  a  smaller  vertical  dimension. 

The  subject  of  diffraction  beam  sidelobe  effects  is  the 
last  remaining  topic  of  this  section.  D.  L.  Hecht  /~22  7  has 
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shown  a  way  to  anticipate  sidelobe  levels,  but  no  one,  to  this 
author's  knowledge,  has  specifically  addressed  sidelobe  sup¬ 
pression.  The  obvious  approach  of  increasing  the  truncation 
ratio  to  suppress  sidelobes  has  the  previously  mentioned 
drawback  of  reduced  information  capacity.  In  present  systems 
having  only  30  dB  overall  dynamic  range,  a  -20  dB  maximum  side¬ 
lobe  level  may  not  be  significantly  harmful.  However,  wider 
dynamic  ranges  in  newer  systems  are  now  evident,  so  eventually 
the  problem  must  be  confronted. 

Since  spatial  sidelobes  are  a  direct  result  of  the  abrupt 
windowing  action  by  the  aperture,  the  only  way  to  remove  them 
is  by  making  the  absorption-transmission  changeover  a  less 
discontinuous  process.  To  produce  a  smoother  window  function, 
a  graded  film  mask  could  be  deposited  on  the  deflector  aper¬ 
ture  exit.  Or,  absorptive  impurities  could  be  diffused  into 
the  interaction  medium  to  produce  the  same  effect.  The  ab¬ 
sorption  profile  would  only  need  to  be  a  function  of  the  hori¬ 
zontal  axis  as  vertical  sidelobes  do  not  interfere  normally  in 
a  one  dimensional  array.  Figure  10  depicts  a  possible  graded 
transmission  profile  which  exhibits  a  smooth  transition  from 
infinite  to  zero  attenuation. 

This  completes  the  formal  discussion  on  laser  beam  pro¬ 
pagation.  To  recapitulate,  a  brief  review  of  propagation 
characteristics  was  presented.  Behavior-determining  equations 
explained  in  the  review  were  utilized  in  a  calculator  program 
for  beam  spot  size  and  radius  prediction.  Gaussian  beam  model 
results  were  compared  with  experimentally  observed  figures  and 
those  predicted  by  a  plane  wave-aperture  interaction  model. 
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D  =  Aperture  Width 


Truncation  ratio  was  found  to  play  a  dominant  role.  Finally, 
the  significance  of  sidelobe  interference  was  discussed  and 
two  methods  of  sidelobe  suppression  were  suggested. 

C.  BRAGG  CELL  PERFORMANCE 

As  mentioned  in  the  introduction,  a  Bragg  cell  translates 
a  time  domain  signal  into  the  frequency  domain  by  converting 
a  signal's  component  frequencies  into  beam  deflection  angles, 
and  respective  intensities  into  beam  intensities.  This  process 
is  sketched  in  Figure  11.  How  well  the  Bragg  cell  accomplishes 
this  task  is  the  subject  of  this  section.  The  three  major 
performance  criteria  of  diffraction  efficiency,  information 
capacity  and  response  time  were  used  to  partition  the  investi¬ 
gation  effort.  The  first  topic  under  analysis  was  diffraction 
efficiency . 

Before  beginning  this  discussion,  a  brief  mention  of  the 
basic  law  governing  Bragg  deflection  is  required.  This  is 
given  by  Equation  (9)  which  relates  the  deflection  angle  to 
the  acoustic  wavelength  of  the  input  signal. 

X 

— -  =  Sin  9=9  (9) 

2A 

o 


or 


whe  re : 

1  =  laser  wavelength 

Aq  =  acoustic  wavelength 


Transform 


focal  focal 

length  length 


n  =  Net  Efficiency 

1  v 


nl  “  n2 


Figure  11.  Bragg  cell  Operation 


9  =  deflection  angle 
f  =  input  frequency 
V  =  acoustic  velocity 


A  detailed  geometric  proof  of  Equation  (9)  is  given  by  Yariv 
in  Ref.  13.  Since  9  is  very  small,  typically  on  the  order 
of  seven  milliradians ,  the  first  order  approximation  is 
sufficiently  accurate  for  spectrum  analyzer  and  optical  ex- 
cisor  applications.  Only  the  two  factors  of  laser  wavelength 
and  acoustic  velocity  are  able  to  alter  the  ratio  of  deflec¬ 
tion  angle  to  input  frequency.  The  laser  wavelength  is  chosen 
from  those  available  to  best  match  the  transmittance  of  the 
Bragg  cell,  and  the  acoustic  velocity  is  a  function  of  other 
constraints,  as  will  be  shown  later. 

On  the  other  hand,  the  beam  diffraction  efficiency  has 
many  factors  which  affect  its  value.  This  parameter  is  de¬ 
fined  as  the  ratio  of  the  diffracted  beam  intensity  to  the 
intensity  of  the  original  undeflected  beam. 

This  is  shown  in  Equation  (10)  which  determines  the  ratio. 
It  was  derived  by  I.  C.  Chang  in  Ref.  21. 


I1  sin2  (n  +  {AI^L/2}2}  ** 

Io  °  (n  +  (AK1L/2}2} 

with 


n 


6  2  PL 

n  P  i  a 
3 '  2 

pV  Hcos  9 
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where: 


1^  =  diffracted  beam  intensity 

Iq  =  original  undeflected  beam  intensity 
=  free  space  laser  wavelength 
n  =  refractive  index 
p  =  elastooptic  coefficient 
p  =  mass  density 
V  =  acoustic  wave  velocity 
=  acoustic  power 
L  =  interaction  length 
H  =  acoustic  beam  height 
AK^  =  wavevector  mismatch 
9q  =  deflection  angle  determined  from  equation  (9) 

An  approximate  version  easier  to  work  with  is  given  by 
Equation  ( 12) : 

I,  0  ,  sin  iAK.L/2} 

yi  =  sin  - - — 5-  (12) 

o  (AKjL/2) 


Let: 

*  .2,1,, 

A  =  sm  in  > 

sin  {AK.L/2} 

B  =  - — -y 

(AK^L/2 } Z 


Term  A  of  equation  (12)  is  constant  with  respect  to  input 
signal  frequency.  Its  effect  will  be  examined  later.  Term  B 
fixes  the  relative  bandshape  of  a  normal  acoustooptic  beam 
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(13) 


deflector.  By  converting  wavevectors  to  frequencies,  B  may 

be  expressed  as: 

sin2  {tt/2L/L  IFF  -F2}} 

B  =  _  o  m 

{tt/2L/L  (FF  -F2}}2 
o  m 

where : 

F  =  normalized  center  frequency  =  1 

F^  =  input  frequency  normalized  to  the  center  frequency 
value 

2 

L  =  nA  cos9  /A  =  characteristic  interaction  lencth 
o  o  o  o 

This  transformation  is  shown  by  D .  L.  Hecht  in  Ref.  20. 

A  plot  of  this  equation  with  various  ratios  of  L/Lo  is  shown 
if  Figure  12.  The  3  dB  bandwidth  is  given  by  equation  (14) 
and  is  also  from  /  20_/. 

Af  =  1.8f  L  /L  (14) 

o  o 

where: 

fQ  =  center  frequency 

Equations  (10)  through  (13)  describe  Bragg  regime  or 
"thick  grating"  diffraction.  Only  a  first  order  beam  will  be 
deflected  with  a  thick  diffraction  grating.  A  thin  diffrac¬ 
tion  grating  will  yield  many  higher  order  diffraction  beams 
in  addition.  This  type  of  diffraction  is  referred  to  as 
Raman-Nath  diffraction.  It  is  intolerable  for  the  intended 
applications  because  of  interference  caused  by  the  additional 
beams . 

Equation  (14)  implies  that  the  deflector  bandwidth  can  be 
easily  increased  by  shortening  the  interaction  length. 
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Figure  12. 


Diffraction  Efficiency  Bandshape 


Increasing  bandwidth  in  this  manner  is  of  limited  overall 
benefit  as  the  diffraction  shifts  from  the  Bragg  regime  to 
Raman-Nath.  In  general,  normal  Bragg  diffraction  is  restric¬ 
ted  to  one  octave. 

This  leaves  only  one  effective  way  in  Equation  (14)  to 
enlarge  the  Bragg  cell's  operating  bandwidth.  The  input 
signal's  center  frequency,  f  ,  must  be  increased.  The  ability 
to  increase  f  is  determined  by  the  input  signal  transducer. 
Substantial  transducer  analysis  is  included  in  Refs.  10  and  21. 
The  highlights  are  briefly  presented  here. 

The  acoustic  transducer  consists  of  a  thin  piezoelectric 
strip  bonded  to  the  interaction  medium.  Its  effectiveness  in 
launching  an  acoustic  wave  into  the  interaction  medium  is  a 
function  of  coupling  and  impedence  matching  characteristics. 

The  particular  bonding  technique  bounds  the  maximum  frequency 
that  the  transducer  can  effectively  couple  to  the  medium. 

The  degree  of  electrical  and  mechanical  impedence  matching 
determines  the  transfer  efficiency  from  electrical  to  acousti¬ 
cal  signal  energy.  It  is  also  frequency  dependent.  Reference 
2  reports  some  available  deflectors  which  operate  in  the 
three  to  four  Gigahertz  region. 

The  above  discussion  of  tranducer  effects  completes  the 
analysis  of  Equation  (12) 's  term  B.  Attention  is  now  directed 
to  term  A.  It  relates  the  interaction  medium's  material  con¬ 
stants  to  diffraction  efficiency. 
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The  argument  of  term  A  may  be  rewritten  as 


h  it  _ 

n  =  j  ‘  c 


where : 


c  =  f—  n3p  /—  • 
\>  J  H 


cosQ 


-}  <  1 


Since  ri’  is  a  sine  function  argument,  improvement  in  diffrac¬ 
tion  efficiency  is  a  "simple"  matter  of  increasing  C.  It  is 
desirable  to  increase  the  refractive  index  -  n  and  elastooptic 
coefficient  -  p,  and  decrease  the  acoustic  beam  height  H, 
the  acoustic  velocity  -  v  and  the  mass  density  -  p.  Little 
can  be  done  with  XQ,  /2  or  9  .  Also,  the  deleterious  effects 
of  increasing  acoustic  input  power  have  been  noted  in  the 
introduction.  And,  inordinate  increasing  of  L  would  prove 
costly  in  terms  of  bandwidth  by  Equation  (14). 

Several  acoustooptic  figures  of  merit  have  gained  accep¬ 
tance  for  use  in  the  selection  of  suitable  interaction  mediums. 
They  are: 


M1  = 


=  M2(nVz} 


6  2 
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(20) 


M4  = 


8  2tt 

n  p  V 


=  M2 { nV2 } 2 


Gordon  introduced  M1  and  M2  in  Ref.  22.  His  results  were 
further  clarified  by  Chang  in  Ref.  21.  Chang  also  included 
a  discussion  of  and  in  equation  (21)  related  diffrac¬ 

tion  efficiency  to  acoustic  power  for  a  given  bandwidth.  This 
is  yielded  by  Equations  (11) ,  (14)  and  (17) . 


M1 

Af 


( 


3 

X  f  H 
o  o 


(21) 


where : 

9  »  1.8tt/2 


This  same  efficiency-power  ratio  can  be  expressed  in  terms 
of  a  fixed  geometry  by  using  M2  in  Equation  (21) . 


is  the  preferred  figure  of  merit  for  systems  in  which  the 
transducer  height  is  not  limited  by  impedence  or  fabrication 
constraints.  H  would  then  be  limited  only  by  the  laser  beam 
height  entering  the  cell.  For  convenience,  H  is  replaced  by 
V/Af,  which  transforms  Equation  (21)  into  Equation  (23). 


f-T-J 

X  f 
o  o 


(23) 
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where : 


P  =  P  /LH  (25) 

a  a 

In  addition  to  having  high  values  of  the  above  figures  of 
merit,  a  good  interaction  medium  should  exhibit  other  charac¬ 
teristics  as  well.  Low  coefficients  of  optical  and  acoustic 
attenuation  are  desirable.  Optical  attenuation  degrades  over¬ 
all  performance  and  ocntributes  heat  to  the  material.  Acoustic 
attenuation  decreases  information  capacity.  Lastly,  if  the 
material  is  crystalline,  then  consideration  must  be  given  to 
how  easily  a  continuous  crystal  structure  can  be  pulled  from  a 
melt.  For  specific  information  on  medium  performance,  see 
Ref.  10  which  contains  an  extensive  comparison  of  material 
properties . 

Once  the  deflector  material  has  been  chosen  and  the  device 
dimensions  fixed,  nothing  remains  free  in  Equation  (10)  to 
improve  diffraction  efficiency  except  for  the  acoustic  input 
power.  However,  beam  steering  and  other  compensation  techni¬ 
ques  can  still  be  called  upon  to  improve  performance  more. 
Korpel's  work  /  5_7  on  stepped  transducer  beam  steering  has 
already  been  briefly  cited.  Gordon  /  22_/  also  discussed  beam 
steering,  only  with  a  planar  instead  of  stepped  array.  Both 
transducers  represented  first  order  steering  as  each  successive 
element  differed  in  phase  by  it  radians. 

Pinnow  /  23_/  compared  normal  or  single  transducer  element 
Bragg  diffraction  with  stepped  and  planar  first  order  beam 
steering.  He  argued  effectively  that  beam  steering  was  better 
as  it  required  less  total  acoustic  power  and  a  lower  cell 
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acoustic  power  density.  This  is  shown  most  clearly  by  attemp¬ 
ting  to  double  bandwidth  in  a  single  element  deflector  cell 
by  doubling  center  frequency.  The  effects  are  revealed  in 
Equation  (14)  with  Lq  rewritten  in  terms  of  frequency. 


L 
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nA  cos0  /A 
o  o  o 
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1 . 8f  L  /L  =  1.8f  K { 
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(27) 


Doubling  the  bandwidth  by  doubling  f  causes  the  transducer 
length  to  reduce  to  one  quarter  as  shown  in  Equation  (28) . 
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This  occurrence  induces  the  acoustic  power  from  Equation  (22) 
to  quadruple. 


~  V 


L/4 

H 


(29) 


The  net  result  is  four  times  the  original  acoustic  power  is 
spread  over  one  fourth  of  the  original  surface  area.  Pinnow 
explained  how  this  unhappy  conclusion  could  be  largely  cir¬ 
cumvented  by  using  a  first  order  array  of  only  moderate  com¬ 
plexity. 

To  begin,  he  pointed  out  that  a  stepped  transducer  array 
is  nearly  twice  as  efficient  as  a  planar  array.  A  first  order 
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planar  array  must  divide  its  power  between  two  beams  which 
are  symmetric  about  the  normal.  The  corresponding  stepped 
array  utilizes  its  inherent  asymmetry  to  put  a  large  percen¬ 
tage  of  the  acoustic  power  into  a  single  main  lobe.  A  typical 
first  order  stepped  array  is  depicted  in  Figure  13. 

The  performance  improvement  provided  by  a  stepped  array 
was  calculated  in  terms  of  an  effective  acoustic  power  ratio. 
To  accomplish  this,  the  array  beam  pattern  needed  to  be 
determined.  For  simplicity,  it  was  described  in  terms  of  far 
field  approximations  even  though  the  light-sound  interaction 
actually  takes  place  in  the  near  field.  The  overall  pattern 
was  taken  to  be  the  produce  of  the  single  element  and  array 
functions  . 


(18)  W(9)  =  B  (9  )  C  (9  ) 


(30) 


where  B(9) 
given  by 

B  (9 )  = 


is  defined  as  the  single  element  function  and  is 


.  - XL9 1 

sini  } 


r  A  L9  , 

I—7-W 


(31) 


C( 9)  is  the  array  function  and  is 


C  (9 ) 


2 

sm 


(9-9')} 


!L£(9-9M  2 


(32) 


9 '  =  beam  steering  angle 
9  =  angle  from  normal  to  element  plane 

D  =  center  to  center  element  spacing 
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Figure  13.  First  Order  Stepped  Array- 
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L  =  element  length 
N  =  number  of  transducer  elements 

The  array  performance  was  compared  at  midband  with  a 
deflector  having  one  element  of  length  Lq  by  taking  the  ratio 
a  of  the  effective  array  beam  power  over  the  effective  single 
element  transducer  beam  power. 


a 


.  LW 
b 

P  L  W 


a  o  s 


(33) 


where : 

I 

=  Array  Acoustic  Power  per  unit  lingth 

=  p'/N 
a 

I 

P  =  Unit  Transducer  Beam  Power  per  unit  length 

cl 

=  Array  Pattern  value  at  midband 
Ws  =  Single  Element  Pattern  at  midband  -  1 


This  leaves 


a 


W  L 
a 

N  L 


W 


a 


3 


(34) 


where : 

3  =  NL/L 

o 

is  the  ratio  of  the  array  length  to  the  single  element  trans¬ 
ducer  lengths. 

At  midband : 

B (9 )  =  G  <  1  (35) 
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2 


2  sin {~{0-9  ' } } 

c  (9  )  =  N  { - - }  (36) 

sin{iD-{0-0  '}} 

»  N2 

therefore : 

a  =  GN  (37) 

The  end  result  is  a  performance  factor  a  equal  to  something 
less  than  the  number  of  transducer  elements.  Pinnow  did  caution 
however,  that  increasing  N  much  beyond  four  was  self-defeating 
as  the  inherent  errors  of  first  order  steering  inhibit  further 
improvement.  Thus,  first  order  beam  steering  can  be  utilized 
to  provide  some  bandwidth  improvement  without  increasing  power. 
Its  benefits  must  be  subjectively  weighed  against  the  diffi¬ 
culties  caused  by  increased  device  complexity. 

Although  beam  steering  does  somewhat  improve  diffraction 
ratio  consistency  across  the  passband,  there  are  also  other 
optical  and  electrical  techniques  which  can  further  minimize 
nonuniformities.  A  set  of  staggered  bandpass  filters  inserted 
prior  to  the  deflector  could  pre-distort  the  xnput  signal  to 
even  out  the  humped  midband  response  of  a  normal  Bragg  cell. 

This  method  would  be  amenable  to  both  spectrum  analysis  and 
optical  excision  systems.  If  advantageous,  an  optical  excisor 
could  also  employ  post-distortion  filtering  equally  as  well, 
however,  a  spectrum  analyzer  would  required  a  scheme  of  time 
varied  gain  for  post-distortion.  An  optical  filter  or  atten¬ 
uator  whose  transmission  is  inversely  related  to  the  diffraction 
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ratio  bandshape  would  be  another  possible  contender.  One 
method  of  obtaining  such  a  filter  employs  photographic  film 
as  the  filter  medium.  The  technique  tailors  the  component  to 
a  particular  system  by  sweeping  a  single  rf  signal  across  the 
deflector  passband  while  the  resulting  diffraction  beam  ex¬ 
poses  an  undeveloped  film  slide  positioned  between  the  de¬ 
flector  cell  and  photodetector.  This  film,  once  developed, 
would  be  a  post-distortion  filter  and  provide  a  uniform  pass- 
band  response,  as  portions  of  the  film  which  were  exposed 
to  higher  illumination  intensities  prior  to  development  would 
attenuate  the  diffraction  beam  to  a  greater  extent. 

The  analysis  up  to  this  point  has  assumed  the  deflector 
treats  each  rf  signal  frequency  component  independently  of 
all  others.  This  supposition  is  not  entirely  true  as  many 
secondary  effects  may  be  noticeably  present.  Each  element 
of  the  composite  input  signal  depletes  the  incident  light 
beam's  capacity  to  respond  linearly  to  the  other  simultaneously 
present  components.  Diffracted  beams  may  also  be  rediffracted 
by  other  acoustic  wavefronts  encountered  prior  to  exiting  the 
cell.  D.  L.  Hecht  /  24_7  analyzed  multiple  frequency  component 
effects.  In  order  to  do  this,  he  used  coupled  mode  equations 
to  solve  the  classical  wave  equation,  Equation  (38)  ,  for  N 
signals  in  both  the  Raman-Nath  and  Bragg  regimes. 


{ u  ( x ,  t )  } 


(38) 


His  solution  algorithm  is  similar  to  that  followed  by  Klein  and 
Cook  in  Ref.  25.  They,  however,  only  treated  the  interaction 
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of  a  single  frequency  sound  column  with  the  incident  light 
beam.  Hecht's  development  and  results  for  Bragg  diffraction 
are  discussed  here  as  this  type  of  interaction  is  more  impor¬ 
tant  for  the  intended  applications. 

The  electric  field  E  was  represented  as  an  N-dimensional 
Fourier  series  expansion  which  was  much  simplified  for  Bragg 
regime  diffraction  since  only  two  principle  modes  were  present. 
This  allowed  the  single  second  order  differential  equation 
expressed  in  Equation  (38)  to  be  broken  up  into  two  first 
order  coupled  mode  differential  equations. 
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dz 
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d  i1  #-* 

Sz  *  (n)  = 


N  V 
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(40) 


where : 

u  =  quiescent  index  of  refraction 

■p°  =  zero  order  diffraction  mode 

^  =  first  order  diffraction  mode 

n  =  (n, ,n0, . n.7)  =  dimensions  of  the  N-tuole  Fourier 

±  £  w 

expansion  of  E  where  n^,n2»  etc. 
have  integer  values  between  -<»  and 


+  °°. 

(n-am)=  <n1'n2' * •nm-i'nm*1'nm+l'nM-l'nN) 
(n+am)=  (nx,n2,  .  •  nn_]_  *  nm+1  *  nm+i '  'nN} 
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Initial  Conditions  at  z  =  0  : 
t°(o)  =  1 
^X(n)  =  0 

Hecht's  results  to  these  equations  are  given  in  terms  of  V 
which  is  "the  normalized  index  of  refraction  amplitude"  for 
the  mth  signal, 
where : 

V  =  Ku  L/cos0  (41) 

m  m  o 


um  =  refractive  index  perturbation  due  to  mechanical  strain 
caused  by  the  mth  signal. 

K  =  2tt/A 

o 


Vm  is  related  to  the  diffraction  ratio  efficiency  quantity 
by  the  following  brief  development. 
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(42) 


where : 


21 


s  = 
m 


m 


PV' 


=  mechanical  strain  due  to  mth  signal 


(Note  that  V  is  the  sound  velocity) 


1^  =  Acoustic  Intensity  due  to  mth  signal 


=  P  /LH 
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Then 
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Equations  (11)  and  (43)  allow  us  to  equate  optical  diffraction 

V 

efficiency  nm  with  -y.  Differentiating  and  decoupling  equations 
(39)  and  (40)  enabled  them  to  be  solved  with  assumed  power 
series  solutions.  The  results  for  N  signals  were  given  in 
terms  of  truncated  series  involving  V_^. 

The  diffraction  efficiency  of  a  single  frequency  component 
signal  in  the  presence  of  IJ-1  other  single  frequency  component 
signals  was  found  to  be  the  following  when  evaluated  at  z  =  L. 
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(44) 


The  corresponding  fractional  deviation  from  ideal  diffraction 
ratio  of  one  signal  due  to  N  signals  has  been  defined  as 
'compression'  C. 


V  2 

T  1  1  T  W  2  M  „  2 
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'Cross  modulation' ,  m,  is  the  compression  of  one  signal  due  to 
a  second  signal  only. 
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2 
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2 


(46) 


(Self  compression  is  not  included  here.) 
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Compression  and  cross  modulation  are  important  considerations 
in  high  density  signal  environments.  Cross  modulation  is  par¬ 
ticularly  so  for  optical  excision  applications,  as  an  amplitude 
modulated  interference  signal  of  high  power  could  do  much  harm 
to  the  reconstructed  signal  of  interest. 

Intermodulation  is  another  area  of  concern  which  is  described 
in  terms  of  a  coupled  mode  solution  ^° ( a a^ ) -  Second  order 
intermodulation  products  are  the  highest  power  intermods  and 
contribute  to  the  zero  order  diffraction  mode.  As  long  as 
I f--f.  <f  .  ,  then  the  second  order  intermod  products  will  not 

be  of  concern.  This  is  guaranteed  if  2f  .  >f 

3  mm  max 

where : 


f  in  =  minimum  deflector  passband  frequency 

f  =  maximum  deflector  cassband  frequency 

max  - 

This  condition  speaks  well  for  deflectors  having  a  bandwidth 
of  slightly  less  than  an  octave.  However,  two- tone  third- 
order  intermods  will  still  fall  within  the  deflector  pass- 
band  and  produce  first  order  diffraction  beams  whose  diffraction 
efficiency  is  described  in  terms  of  -p  (Za^-a^)  as  follows: 


1  1  V1  4  V2  2 
X2,-1  =  *  <2ar*2)  =  IT  {-r}  {~r] 


(47: 


Three-tone  third-order  intermods  of  the  form  (f^+f^-f^)  will 
also  create  first  order  diffraction  beams  in  accordance  with: 


I 


1,1, 


-1 


(al+a2"a3) 


(48) 
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As  system  dynamic  ranges  improve,  third  order  intermodulation 
products  will  continue  to  become  more  significant. 

These  last  remarks  complete  discussion  on  multiple  fre¬ 
quency  performance  effects.  Attention  is  now  drawn  to  another 
closely  related  subject  concerning  geometric  implications  of 
light  and  sound  interaction.  Uchida  /  10_/,  in  his  derivation 
of  diffraction  efficiency,  assumed  weak  interaction  of  plane 
wave  light  incident  to  a  rectangular  sound  column.  The  re¬ 
sults,  while  very  useful,  do  not  take  into  account  true  geo¬ 
metric  considerations  of  the  light  and  sound  beams.  Several 
authors  have  challenged  various  aspects  of  the  problem.  Two 
particularly  illuminating  papers  are  briefly  discussed  here. 
The  first  paper  /  26_/  presented  closed-form  or  analytic  solu¬ 
tions  to  problems  where  the  light  and  sound  beams  were  of 
similar  profiles. 

The  most  common  type  of  light-sound  interaction  experi¬ 
mentally  observed  is  between  a  Gaussian  profile  light  beam 
and  a  rectangular  sound  column.  McMahon  /  26_7  argued  that 
a  closed-form  description  of  this  instance  could  not  be 
arrived  at.  Instead,  he  proposed  to  bracket  anl  indirectly 
analyze  this  case  by  finding  solutions  for  two  closely  related 
cases.  The  two  he  chose  were  a  light  beam  of  rectangular 
profile  intersecting  a  like  sound  column  and  a  Gaussian  to 
Gaussian  intersection.  He  conducted  the  analysis  by  using 
an  oscillating  electric  dipole  model  formulation  and  was  able 
to  determine  the  first  order  beam  power  for  both  situations. 
The  following  results  were  derived  for  the  light  beam  incident 
at  the  Bragg  angle  so  that  wavevector  mismatch  would  not  have 


to  be  considered. 


For  the  rectangular  to  rectangular  interaction 


rr 


=  A{  1 


-!> 


(49) 


where 


P  =  power  in  first  order  diffraction  beam  for  rectan- 
rr 

gular  to  rectangular  interaction 


A 


=  f(Pj 

a. 


2 

WL 

cos  9 

o 


f(P  )=  function  of  input  acoustic  power  and  material 

cl 

constants 

W  =  laser  beamwidth 

L  =  sound  column  width 

2Lsin0 

a  =  -rr  o 


The  Gaussian-Gaussian  interaction  gives 

Pgg  =-it2A/{1  +  a2}*5  (50) 

where : 

P  =  power  in  first  order  diffraction  beam  for  Gaussian 

gg 

to  Gaussian  interaction 

The  author  claimed  for  the  condition  of  (Jj<a<2)  that  the 
power  diffracted  by  either  configuration  would  vary  from  the 
other  by  no  more  than  eight  percent.  Thus,  by  taking  the  mean 
value  between  them,  one  could  only  be  off  by  about  four  percent 
from  the  "true"  value  for  the  intermediate  Gaussian  to 
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rectangular  interaction  case.  However,  for  the  deflector  used 
in  the  experimental  set  up  of  this  thesis,  the  calculated 
figures  for  the  two  configurations  differed  by  a  considerable 
margin  since  a  was  not  within  the  specified  limits. 

W  ~  L 

0  =7  milliradians 

o 

a  =  .014 

P  =  ( .  995) A 

rr 

P  =  {9 . 86) A 

gg 

In  spite  of  this  variation,  the  results  do  give  some  insight 
to  geometric  influences  on  diffraction  efficiency.  To  obtain 
a  more  general  appraisal  of  the  interaction  process,  Korpel's 
work  in  Ref.  27  should  be  consulted.  He  presents  a  technique 
for  handling  light  and  sound  beams  of  arbitrary  profile.  The 
ensuing  derivation  assumed  a  spatially  varying  plane  wave 
angular  spectrum  representation  of  the  radiated  fields  and  a 
"multiple  scattering  model"  for  the  interaction  process. 

Numerical  results  for  particular  beam  configurations  require 
power  series  solutions  that  the  author  left  for  future  deter¬ 
mination. 

In  conclusion,  primarily  qualitative  benefits  are  available 
from  McMahon's  analysis  of  geometric  factors.  And  while  Korpel's 
theory  offers  exact  solutions  for  general  beam  profiles,  im¬ 
plementation  may  be  difficult  in  specific  cases.  Much  further 
work  remains  in  this  area. 

The  next  and  final  sub-topic  related  to  diffraction  effi¬ 
ciency  is  a  discussion  of  Bragg  cell  dynamic  range.  This 
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performance  measure  is  defined  as  the  linear  operating  range 
limited  by  optical  background  noise  for  minimum  detectable 
signals  and  saturation  for  large  signals.  There  are  three 
types  of  optical  scattering  which  are  causes  for  background 
noise.  Airborne  dust  and  lens  defects  are  scattering  centers 
included  in  the  first  class.  Their  noise  contribution  is 
basically  a  function  of  the  environment  and  component  quality 
for  a  particular  operating  system,  and  so  will  not  be  consi¬ 
dered  here.  By  comparison,  optical  scattering  resulting  from 
the  diffraction  of  the  zero  order  beam  is  a  function  of  design 
and  is  one  of  the  two  possible  theoretical  limits  for  minimum 
background  noise.  The  other  possible  limit  is  scattering 
due  to  thermally  generated  acoustic  waves  in  the  Bragg  cell 
medium.  This  type  of  scattering  in  known  as  Brillouin 
scattering.  The  saturation  limit  for  large  signals  is  deter¬ 
mined  by  zero  order  depletion  considerations  in  multiple 
frequency  environments.  Another  maximum  linear  response  limit 
could  be  used  in  signal  environments  composed  of  only  one 
frequency  component  at  a  time,  however  this  limit  would  have 
no  practical  value  for  typical  operating  conditions. 

To  determine  a  system's  minimum  detectable  signal,  both 
theoretical  limits  must  be  calculated  and  compared  for  the 
larger  value.  The  minimum  detectable  signal  due  to  zero  order 
diffraction  effects  is  considered  first.  In  the  worst  case 
(maximum  sidelobes) ,  the  diffraction  pattern  of  the  main  beam 
obeys  Equation  (7) .  Sidelobe  maxima  are  roughly  determined 
by  /  20  7’ 
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(51) 
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where : 

x  =  D0/X 

o 

N  =  angle  of  deflection  due  to  the  minimum  passband  fre¬ 
quency/angular  width  of  a  single  spot  and  is  given 
by  D9min/Ao. 

The  point  where  the  first  order  beam  corresponding  to  the 
minimum  passband  frequency  intersects  the  zero  order  dif¬ 
fraction  pattern  determines  a  minimum  detectable  signal  level 
due  to  zero  order  interference  effects.  This  minimum  signal 
value  relative  to  the  axis  intensity  of  the  zero  order  beam 
was  calculated  as  follows  for  the  NPS  a-o  spectrum  analyzer 
assuming  an  octave  bandwidth: 

Let 


I  =  I  (9  .  )  =  minimum  detectable  signal 
s  o  min 

Then 
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D  =  7  mm 


72 


It  is  important  to  note  that  the  above  calculated  figure  is 
independent  of  wavelength.  The  only  way  to  lower  minimum 
detectable  signal  level  due  to  diffraction  is  by  sidelobe 
reduction.  For  systems  in  which  sidelobes  are  sufficiently 
attenuated,  the  minimum  required  signal  would  be  limited  by 
Brillouin  scattering.  The  scattered  power  due  to  this  effect 
is  given  in  the  following  equation  /~28_7: 


^T  =  8ttFL 

P  3 
o 

where : 


(53) 


PT  =  total  power  scattered  in  all  directions 
Pq  =  incident  power 

F  =  scattering  coefficient  proportional  to  T3 
T  =  temperature  in  degrees  Kelvin 
L  =  interaction  medium  length 
In  the  subsequent  calculation  of  minimum  ratio,  F  was  not 
available  for  Lead  Molybdate  so  a  representative  coefficient 
for  a  very  strong  scatterer  was  used. 


2FLS 

j  =  normalized  part  of  incident  power  delivered 


into  one  resolution  cell 
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3  -130  dB 

where: 

F  =  40  cm"1 

L  =  .7  cm 

S  =  A0 (1  -  cosA 0 ) 

A0  =  \q/L 

\  =  .6328  um 

o 

If  the  two  power  and  intensity  ratios  are  taken  to  be  essen¬ 
tially  equivalent  measures,  then  the  zero  order  diffraction 
sidelobe  level  was  clearly  much  stronger  than  the  noise  level 
due  to  Brillouin  scattering.  The  upper  linearity  limit  was 
determined  by  the  maximum  single  frequency  input  allowable 
in  a  multifrequency  environment  which  has  been  shown  by  Hecht 
/  24 _7  to  be: 

ni<  .01  =  -20  dB 

The  resulting  dynamic  range  (DR)  calculated  for  the  NPS  Lab 
deflector  with  this  condition  was  then: 

DR  =  24  dB 

In  general,  it  would  appear  from  the  above  typical  calcu¬ 
lations  that  systems,  which  do  not  employ  special  sidelobe 
attenuation  techniques,  would  be  zero  order  diffraction 
limited  in  dynamic  range.  However,  F  is  very  temperature  de¬ 
pendent  so  a  unit  might  begin  operation  limited  by  zero  order 
sidelobes  and  gradually  suffer  a  reduction  in  dynamic  range 
due  to  device  heating. 


With  discussion  of  diffraction  efficiency  now  completed, 
concentration  is  shifted  to  information  capacity.  This  per¬ 
formance  measure  is  defined  as  the  number  of  beam  spots 
resolvable  in  the  image  plane  and  is  determined  by: 


20  =  Xq/A  =  deflection  angle  from  zero  order  beam 


A4>  =  2(9 


max 


=  X  Af/V 
o 

N  =  A<J>/A0 


9  .  ) 
mm 


=  AfT  =  number  of  resolvable  beam  spots  in  image  plane 

where : 

x  =  D/V 

D  =  effective  aperture  width 
V  =  acoustic  velocity 

A  large  information  capacity  is  desirable  but  difficult  to 
obtain,  as  only  two  factors  affect  its  value.  Bandwidth 
improvement  considerations  have  already  been  discussed  at 
length  under  diffraction  efficiency  and  so  are  not  covered 
here.  The  length  of  time  x  required  for  a  wavefront  to  travel 
the  effective  length  of  the  sound  column  is  the  principal 
topic  of  this  section.  To  lengthen  x,  a  low  acoustic  velocity 
would  apparently  assist  matters,  but  acoustooptic  materials 
which  have  low  sound  velocities  also  tend  to  have  high  acoustic 
loss  factors.  This  attenuation  limits  the  deflector's  effec¬ 
tive  aperture  width  and  gives  the  sound  column  an  exponential 
amplitude  profile  shown  by: 


A  (x) 


g-ax/V 


(55) 
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where : 


A (x)  =  normalized  acoustic  amplitude 

x  =  distance  from  transducer  along  sound  column 

a  =  acoustic  loss  factor 

2 

=  Tf  V/{201og^Qe}  nepers/sec 

2 

T  =  material  acoustic  loss  constant  (dB/{cm-GHz  .r ) 

V  =  acoustic  velocity 

In  order  to  gain  a  rough  measure  of  the  effective  aperture's 
relation  to  acoustic  loss  factor,  an  arbitrary  width  was 
chosen  where  A(D)  equals  2/e  or  approximately  1//2.  Brief 
algebraic  manipulation  yields: 

D  =  {Vloge2}/a  (56) 

=  6/{Tf2} 


The  number  of  resolvable  frequencies  may  now  be  rewritten  as: 


N  =  Kf  /f2;  let  f  =  f 
o  o 

=  K/fo 


(57) 


where : 

K  =  llLQ/{vrL} 


This  result  clearly  shows  better  acoustic  materials  have 
smaller  VT  products.  It  also  points  out  for  a  given  medium, 
maximum  information  capacity  is  inversely  proportional  to 
center  frequency.  The  implication  for  operating  systems  is 
that  once  the  "best"  material  has  been  selected  for  a  normal 
diffraction  Bragg  cell  then  only  a  reciprocal  relationship 
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exists  between  the  system's  surveillance  bandwidth  and  infor¬ 
mation  capacity.  One  must  be  sacrificed  by  the  same  percentage 
the  other  is  improved.  However,  this  last  conclusion  repre¬ 
sents  the  ideal  case,  inferior  performance  would  be  observed 
for  systems  which  allow  thermal  gradients  to  form  during 
periods  of  high  power  input.  In  Ref.  10,  Uchida  presented 
findings  on  this  subject  which  showed  in  one  case  as  much  as 
a  five  percent  drop  in  information  capacity  for  an  input  power 
increase  of  six  percent.  Proper  heat  sinking  was  reported 
to  significantly  reduce  this  heat  dependence. 

Quick  response  is  the  last  general  topic  in  the  discussion 
of  Bragg  cell  performance  factors.  This  performance  measure 
is  desirable  in  some  applications,  such  as  Bragg  cells  intended 
to  be  beam  modulators  rather  than  deflectors.  On  the  other 
hand,  high  resolution  surveillance  spectrum  analyzers  would 
not  be  noted  for  response  agility  as  information  capacity  would 
directly  suffer.  A  txadeoff  must  be  made  for  optical  excisors 
as  both  resolution  and  response  speed  are  desirable  here.  For 
the  purpose  of  this  work,  the  sharpness  of  Bragg  cell  response 
is  defined  as  one  over  t  — the  rate  at  which  single  resolvable 
frequency  inputs  may  be  "hopped"  and  still  attain  maximum  in¬ 
tensity  diffraction  ratio.  It  is  a  measure  of  how  "instan¬ 
taneous"  the  displayed  spectrum  is.  The  diffraction  intensity 
of  short  pulses  is  determined  by  /— 29 _7: 

D 

Is(t)  =  /oA(x,t)IQ(x)dx  (58) 
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where: 


Is(t)  =  "Short  Pulse"  instantaneous  light  intensity 
response 

A (x , t)  =  relative  diffraction  of  signal  in  medium  as  a 
function  of  transited  distance  and  time 
I  (x)  =  incident  light  beam  cross  section  in  aperture 

D  =  effective  aperture  width 

Short  pulse  attenuation  effects  were  particularly  noticeable 
when  observing  the  spectrum  of  a  sinusoidally  modulated  FM 
imput  signal  shown  in  Figure  14.  The  input  power  level  re¬ 
quired  an  increase  of  nearly  two  orders  of  magnitude  in 
contrast  to  the  power  needed  for  a  single  beam.  For  a  quali¬ 
tative  overall  performance  comparison,  a  display  of  the  same 
signal  spectrum  produced  by  a  conventional  scanning  super¬ 
heterodyne  spectrum  analyzer  is  presented  in  Figure  15. 

In  this  section,  a  thorough  compilation  and  examination 
of  all  important  parameters  was  attempted.  To  summarize  the 
findings:  large  information  bandwidth,  quick  response,  high 

diffraction  efficiency  and  large  dynamic  range  are  functions 
of  geometry,  material  constants  and  tradeoffs. 

D.  PHOTODETECTOR  PERFORMANCE 

The  Fairchild  CCD121H  and  Reticon  CCPD-1728  photodetector 
arrays  were  the  principal  detectors  examined  in  this  work  for 
spectrum  analysis  application.  The  two  1728  photosite  arrays 
have  many  similar  features  of  interest,  especially  in  their 
internal  working  structures.  Each  detector  photosite  of  both 
arrays  integrates  the  current  generated  by  the  incident  light. 
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Figure  14.  NPS  A-0  Spectrum  Analyzer 

Display  of  Wideband  FM  Signal 


/ 

i  I- 


Figure  15.  Conventional  Spectrum  Analyzer 
Display  of  Same  Wideband  FM 
Signal  as  in  Figure  14 


For  this  charge  storage,  the  Reticon  array  photosites  use 
photodiode  junction  capacitance  while  the  Fairchild  device 
elements  rely  on  dielectric  capacitance.  Both  arrays  transfer 
the  stored  charge  to  CCD  analog  shift  registers  for  further 
transfer  to  amplifying  circuits  and  eventual  display.  This  is 
accomplished  by  dividing  the  array  elements  into  even  and  o id 
groupings.  Even  elements  transfer  to  a  shift  register  on  one 
side  of  the  photosite  array  and  odd  ones  transfer  to  an 
opposite  register.  The  first  difference  in  operation  occurs 
when  the  registers  are  transferring  data  off  chip.  The  CCD 
chip  interleaves  the  even  and  odd  charge  packets  prior  to  off 
chip  exit.  The  CCPD  ship  does  not,  but  rather  leaves  this 
task  for  external  circuitry.  Diagrams  of  both  chip  operations 
are  shown  in  Figures  16  and  17. 

Examination  of  respective  data  sheets  for  the  two  devices 
reveals  near  equivalence  in  all  specifications  save  two.  The 
CCD  array  was  a  denser  and  faster  chip.  It  required  fourteen 
percent  less  area  per  photosite  and  used  rectangular  vice 
square  shaped  elements  with  the  short  dimension  along  the 
array  face.  This  translated  into  a  twenty  percent  resolution 
improvement  over  the  photodiode  chip.  Also,  the  maximum  output 
data  rate  of  the  Fairchild  device  was  twice  that  of  the  Reticon 
chip.  This  meant  better  resolution  of  input  signal  temporal 
characteristics . 

Dynamic  range  for  both  chips  was  defined  as  the  individual 
array  element's  ratio  of  saturated  output  voltage  to  background 
noise  equivalent  voltage.  The  advertised  dynamic  range  for 
each  device  was  given  as  27  dB.  This  figure  was  verified  by 


Figure  16.  Fairchild  Charge  Coupled  Device 
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Figure  17.  CCPD  Photodetector  Array  Chip  Schematic 


two  separate  methods  for  the  NPS  spectrum  analyzer  CCD 
detector.  The  first  method  utilized  the  calibration  marks 
on  the  rf  input  signal  generator  for  this  measurement.  The 
subsequently  observed  dynamic  range  was  27.5'  dB  plus  or  minus 
one  dB.  The  second  technique  employed  neutral  density 
filters  to  attenuate  the  saturated  output  signal  down  to  the 
noise  level  and  yielded  a  26  dB  range.  However,  slight  beam 
spreading  due  to  filter  presence  was  noted  which  further 
attenuated  the  signal  and  accounted  for  the  lower  assessment. 

Brief  experience  was  gained  with  a  smaller  Reticon  CCPD- 
256  chip  and  related  support  circuitry  before  working  with 
the  Fairchild  CCD-1728  chip  and  development  board.  The  photo¬ 
diode  array  required  slightly  less  effort  to  fine  tune  than 
the  Fairchild  model,  however  this  may  have  been  due  to  the 
simpler  nature  of  the  smaller  device  package.  The  larger 
Fairchild  unit  was  used  extensively  over  an  eight  month  period 
and  performed  with  100  percent  reliability  throughout  this 
time.  The  author  believes  either  device  technology  would 
perform  commendably  in  spectrum  analyzer  applications,  but  is 
biased  in  favor  of  the  straight  CCD  approach  in  view  of  the 
above  noted  advantages. 

Optical  excisors,  as  noted  in  the  introduction's  list  of 
specific  goals,  require  a  detector  much  different  from  the 
devices  discussed  above.  Most  excisor  applications  require 
wide  information  bandwidths  of  the  kind  offered  by  photomul¬ 
tiplier  tubes  (PMT)  available  on  an  "off  the  shelf"  basis. 

For  example,  Varian  markets  typical  PMTs  which  have  bandwidths 
of  a  GigaHertz  and  dynamic  ranges  of  114  dB .  These  specifications 
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mean  that  a  given  system's  performance  will  be  most  likely 
limited  by  the  Bragg  cell  dynamic  range  and  bandwidth.  With 
these  comments,  focus  is  now  placed  on  performance  modeling 
of  the  optical  excision  process. 

E.  ELECTRONIC  AND  OPTICAL  EXCISOR  PERFORMANCE  MODELS 

Optical  excision  performance  is  measured  in  terms  of  two 
gauges  which  are:  how  well  offending  signal  frequency  com¬ 
ponents  are  removed  from  the  output,  and  how  much  the  original 
signal  of  interest  is  corrupted  by  the  loss  of  frequency  com¬ 
ponents  shared  with  the  offending  signal.  Both  gauges  are 
primarily  functions  of  the  electronic  or  optical  clipper 
selected  for  use  in  the  particular  system.  In  this  analysis, 
electronic  clippers  are  assumed  to  completely  remove  inter¬ 
ference  signals  chosen  to  be  excised  since  nearly  infinite 
attenuation  would  be  available  at  the  output  ports.  This 
same  claim  cannot  be  made  for  optical  clippers  because  of 
diffraction  limitations.  The  narrow  obstacle  of  the  optical 
clipper  impedes  the  progress  of  the  intended  beam  in  much  the 
same  manner  as  a  tree  or  telephone  pole  positioned  between 
two  people  impedes  dialogue.  The  observed  diffraction  effects 
are  advantageous  for  conversation  but  not  for  excision. 

Examination  of  optical  clipper  efficiency  was  approached 
by  assuming  the  obstruction  was  much  narrower  than  the  im¬ 
pinging  diffracted  beam.  In  the  limit,  this  situation  may 
be  described  as  a  narrow  obstacle  in  a  coherent  plane  wave 
light  field.  Clipper  effectiveness  then  becomes  a  function 
of  the  narrow  obstacle  width  and  the  distance  to  the  image 
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plane  depicted  in  Figure  18.  The  resulting  image  plane  in¬ 
tensity  diffraction  pattern  provides  a  visual  measure  of 
clipper  utility.  This  diffraction  pattern  was  calculation  in 
terms  of  Fresnel  integrals  where  the  narrow  obstacle  was 
considered  to  be  the  intersection  of  two  overlapping  infinite 
half  planes.  For  further  information,  brief  developmental 
guidance  to  derive  the  following  image  plane  intensity  profile 
is  given  in  Ref.  17.  The  normalized  light  intensity  profile, 

I (U)  ,  along  the  u  axis  in  the  image  plane  (Fig.  19)  is  given 
by: 

I(u)  =  ( (B  -  A)2  +  (2  -  A  -  B)2}/4  (59) 

where : 

U1  2 

A  =  /  costtw  /2  dw 

u2 

U1  2 

B  =  /  sinirw  /2  dw 

u2 

2  L 

u2  =  ^2  ^X — '  =  normalizec^  integral  limit 

o 

y 2=  actual  integral  limit  mapped  to  u^ 

Ay=  Y i  ~  ¥ 2  =  actual  obstacle  width 

rQ=  distance  from  obstacle  to  image  plane 

X  =  light  wavelength 

This  equation  was  approximated  by  Simpson's  rule  and  used 
in  a  Fortran  program  to  plot  image  plane  intensity  profiles 
shown  in  Figure  20  through  26.  The  figure  sequence  shows  the 
effect  of  an  obstacle  which  is  moved  steadily  away  from  the 
image  plane  until  its  effect  is  negligible.  The  program 
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Figure  18.  Narrow  Obstacle  in 


Plane  Wave  Light  Field 
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in  front  of  the  image  plane. 


for  a  16  um  wide  obstacle  located  10 
in  front  of  the  image  plane. 
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Figure  26.  Image  plane  intensity  profile  vs  distance 
for  a  16  um  wide  obstacle  located  1.2m 
in  front  of  the  image  plane. 


used  for  plotting  is  listed  in  Appendix  C. 

To  analyze  the  effects  of  signal  corruption  by  the  excisor, 
the  reader  is  directed  to  consult  Ref.  30.  In  this  work, 
Sussman  and  Ferrari  conducted  a  study  of  Direct  Sequence 
Spread  Spectrum  correlation  peak  attenuation  due  to  the  pre¬ 
sence  of  an  inband  notch  filter.  They  found  that  interference 
frequency  spikes,  which  have  bandwidths  of  up  to  three  percent 
of  the  spread  spectrum  signal  bandwidth,  can  be  excised  with 
only  one  dB  drop  in  the  output  correlation  peak.  This  worst 
case  result  was  obtained  for  a  notch  filter  centered  on  the 
spread  spectrum  signal  center  frequency.  Much  proportionately 
wider  excision  bandwidths  were  tolerated  for  slightly  offset 
notch  filter  center  frequencies.  In  conclusion,  their  analysis 
indicated  the  comparatively  slight  signal  degradation  due  to 
the  notch  filter  was  much  ^referable  to  the  degradation  caused 
by  the  narrowband  interference. 
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III.  SUMMARY/RECOMMEND ATI ONS 

A.  RESULT  SUMMARY 

Laser  beam  propagation,  Bragg  cell  performance,  photo¬ 
detector  performance,  and  optical  excision  models  were  the 
main  topics  addressed  in  this  work.  For  the  laser,  optimal 
light  beam  profile  behavior  in  the  image  plane  was  found  to 
be  a  function  of  both  Gaussian  beam  propagation  and  aperture 
truncation  effects.  Observations  of  an  assembled  spectrum 
analyzer  were  made  which  compared  closely  with  calculated 
values.  And,  two  methods  of  beam  sidelobe  attenuation  were 
proposed  for  signal  resolution  improvement.  Bragg  cell  per¬ 
formance  was  measured  in  terms  of  diffraction  efficiency, 
information  capacity,  and  response  agility.  Also,  device 
dynamic  range  was  calculated  for  a  multifrequency  input  signal 
environment.  Additionally,  beam  steering  and  diffraction 
efficiency  effects  of  light-sound  beam  geometric  factors  were 
studied.  Photodetector  array  performance  of  charge-coupled  and 
photodiode  technologies  was  compared  for  spectrum  analyzer 
application.  Optical  excisor  models  were  discussed  in  terms 
of  interference  removal  effectiveness  and  subsequent  corruption 
of  residual  information  bearing  signals.  A  Fresnel  approximation 
model  of  optical  clipping  was  developed  and  programmed  for  com¬ 
puter  plot  generation.  And  finally,  signal  degradation  was  dis¬ 
cussed  for  DSSS  signal  correlation  peak  attenuation  due  to 
excision. 
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B.  RECOMMENDATIONS 


The  following  remarks  are  specific  recommendations  for 
optimal  performance  of  spectrum  analysis,  and  optical  excision 
of  narrowband  interference.  To  reduce  zero  order  sidelobe 
interference,  and  thereby  increase  dynamic  range,  use  of  Bragg 
cell  apertures  which  do  not  abruptly  change  in  optical  trans¬ 
mittance  at  the  window  edges  is  recommended.  Choice  of  a 
deflector  interaction  medium  having  high  refractive  index, 
large  elastooptic  coefficient,  small  mass  density,  and  small 
acoustic  velocity-material  acoustic  loss  constant  produce  is 
also  desirable.  Either  photodetector  technology  was  found 
acceptable  with  a  slight  advantage  held  by  the  charge-coupled 
device  array.  Unfortunately,  no  guidance  can  be  given  with 
regard  to  the  choice  of  a  "best"  clipper  for  optical  excision 
as  developments  in  this  area  are  still  progressing  quickly. 

C.  AREAS  FOR  FURTHER  ANALYSIS 

Further  effort  related  to  this  thesis  is  required  in 
Gaussian  beam  truncation  analysis  and  electronic/optical 
clipping  model  development.  A  Discrete  Fourier  Transform 
(DFT)  computer  model  of  aperture  windowing  of  Gaussian  beams 
is  one  possible  way  of  promise  to  study  truncation  effects  as 
analytic  solutions  are  exceptionally  difficult  to  arrive  at. 
Determination  of  window  functions  which  "best"  attenuate 
sidelobes  and  maximize  information  capacity  would  be  highly 
desirable.  Much  further  work  also  remains  in  optical  excision 
modeling.  Another  DFT  computer  simulation  could  be  used  to 
render  a  Fraunhofer  apploximation  model  of  obstacle-light  beam 
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interaction  in  order  to  further  clarify  optical  clipper 
effectiveness.  This  approach  could  easily  handle  the  case 
where  the  obstacle  width  would  be  equal  to  the  incident  light 
beam  width.  Lastly,  a  similar  model  could  be  used  to  determine 
nearly  exact  degradation  effects  of  the  excision  process.  These 
results  could  then  determine  optimal  utilization  of  the  exci¬ 
sion  process. 

D.  CONCLUSION 

The  performance  analysis  guidelines  developed  in  the  body 
of  this  thesis  provide  a  realistic  view  of  acoustooptic  poten¬ 
tial  for  spectrum  analysis  and  interference  excision  applica¬ 
tions.  Also,  several  acoustooptic  spectrum  analyzers  have 
been  made  with  commercially  available  components  and  readily 
provide  wideband,  high  resolution  coverage  not  easily  obtained 
in  conventional  systems.  Lastly,  work  in  optical  excision 
techniques  has  progressed  substantially  and  ongoing  efforts 
continue  to  show  performance  improvement.  In  view  of  these 
facts,  further  research  effort  in  this  area  is  very  likely  to 
provide  many  more  valuable  benefits  for  applications  involving 
wideband,  high  density  signal  environments. 
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APPENDIX  A 


< 


BEAM  PHASE  FRONT  RADIUS  AND  SPOT  SIZE  PLOT  PROGRAM 

This  Fortran  program  plots  Gaussian  beam  phase  front 
radius  and  spot  size  as  a  function  of  distance  from  beam 
waist.  Selectable  parameters  are: 

LMBDA  =  laser  wavelength 

NPTS  =  number  of  plot  points  desired 

WZERO  =  waist  beam  spot  size 

DCRM  =  distance  increment 

Prospective  users  at  the  Naval  Postgraduate  School  should 
refer  to  the  W.  R.  Church  Computer  Center's  Versatec  plotter 
routine  user  manual  for  required  job  control  cards  and  plot 
dimensions . 


f 
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C  BEAM  RADIUS  AND  WAIST  AS  A  FUNCTION  OF  DISTANCE 

DIMENSION  R(500) ,  X(500),  W(500) 

INTEGER* 4  ITB (12) /12*0/ 

REAL*4  RTB (28)/28*0.0/ 

RPTS  =  500 
REAL  LMBDA 
LMBDA  =  .6328  E-06 
WZERO  =  .44  E-03 
LMT  =  NPTS  +  1 
DCRM  =0.03 
PI  =  3.141592654 
1  =  1 

Z(l)  =  1.0  E-03 

100  R ( I)  =  Z (I)  +  (PI  *  WZERO* *2/LMBDA) **2/Z (I) 

A  =  SQTR (1.0  +  (LMBDA  *  Z (I) / (PI*WZERO**2) ) **2) 
W(I)  =  WZERO  *  2 
I  =  I  +  I 

IF  (I  .EQ.  LMT)  GO  TO  999 

Z (I)  =  Z (I  -  1)  +  DCRM 

ITB (3)  =  8 

ITB  (4)  =  4 

RTB ( 1 )  =  .5 

RTB (2)  =  .5  E-03 

GO  TO  100 

999  CALL  DRAWP  (NPTS , Z ,W, ITB , RTB) 

RTB (2)  =5.0 

CALL  DRAWP  (NPTS , Z ,R , ITB , RTB) 

STOP 

END 
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GAUSSIAN  PROPAGATION  PROGRAM 

This  TI-58,  TI-59  hand-calculator  program  calculates 
Gaussian  beam  phase  front  radius  and  spot  size  for  propaga¬ 
tion  through  any  number  of  thin  lenses  spaced  at  arbitrary 
distances.  Selectable  parameters  are  stored  in  the  following 
registers : 

*R00  =  laser  wavelength 
R01  =  beam  spot  size 
R03  =  lens  focal  length 
R04  =  distance  increment 

Intermediate  results  are  stored  in: 

R06  =  real  part  of  complex  q 

RO 7  =  imaginary  part  of  complex  q 

*  Common  units  are  used  throughout  the  program 

OPERATING  SEQUENCE 

The  wavelength  of  the  laser  is  the  first  necessary  entry. 
This  value  is  entered  into  program  memory  by  loading  it  into 
the  display  register  and  pressing  label  A.  The  initial  beam 
spot  size  and  phase  front  radius  are  entered  by  first  keying 
spot  size  into  the  display  and  pressing  B.  Immediately  after 
this,  the  phase  front  radius  is  loaded  in  the  display  and 
entered  by  pressing  the  R/S  key.  If  the  user  desired  to  start 
at  the  waist  then  an  infinite  phase  front  radius  can  be 
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entered  by  loading  any  number  greater  than  or  equal  to  1000, 
as  all  such  values  for  phase  front  radius  are  treated  as 
infinite. 

The  above  portion  of  the  sequence  determines  the  initial 
conditions  for  beam  propagation.  From  this  starting  point, 
any  amount  of  propagation  distance  may  be  entered  by  loading 
the  additional  path  length  desired  in  the  display  and  pressing 
D.  The  newly  calculated  beam  spot  size  is  stored  in  the 
display  register  and  the  corresponding  phase  front  radius  is 
stored  in  the  t-register.  Propagation  through  thin  lenses  is 
handled  in  a  similar  manner  by  entering  the  focal  length  in 
the  display  and  pressing  C.  Spot  size  and  phase  front  radius 
are  read  out  in  the  same  manner  as  for  additional  propagation 
distance. 

The  above  comments  complete  the  operating  sequence  des¬ 
cription.  This  is  now  followed  by  a  listing  of  the  program 
steps  by  their  functional  sections. 

SEQUENCE  FOR  ENTERING  LASER  WAVELENGTH: 

*Lbl  A  STO  00  R/S 

SEQUENCE  FOR  ENTERING  INITIAL  BEAM  SPOT  SIZE  AND  PHASE  FRONT 
RADIUS : 

"Lbl  B  STO  01  R/S  STO  02  1000  x<>t  RCL  02  *x>t  INV 
(RCL  02  1/ x  t  RCL  02  x2  1/x  +  (RCL  00  r  *tt  t  RCL  01 
x2) )  x2) )  STO  05 

(RCL  00  r  *tt  t  RCL  01  X2  t  (RCL  02  1/x  X2 

T  *n  T  RCL  01  X2)  X2))  STO  06 
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+  (RCL  00 


CLR  X<>t  CLR  R/S 

*Lbl  INV  CLR  STO  05  (RCL  01  X2  x  r  RCL  00) 

STO  06  CLR  X<> t  CLR  R/S 

SEQUENCE  FOR  PROPAGATION  THROUGH  SPATIAL  DISTANCE: 

*Lbl  D  +  RCL  05  =  STO  05  E 

SEQUENCE  FOR  CALCULATING  BEAM  SPOT  SIZE  AND  PHASE  FRONT 
RADIUS  FROM  COMPLEX  "q" : 

*Lbl  E  ((RCL  05  x2  +  RCL  06  x2)  r  RCL  05)  STO  02 
xot  (RCL  00  t  *ir  t  RCL  06  x  (RCL  05  x2  +  RCL  06  x2 
)  )  /x  STO  01  R/S 

SEQUENCE  FOR  TRANSMISSION  THROUGH  THIN  LENS: 

*Lbl  C  STO  03  (RCL  03  x2  x  RCL  06  t  ((RCL  03  - 

RCL  05)  x2  +  RCL  06  x2  ))  x<>t  ((RCL  03  X  RCL  05 
x  (RCL  03  x  RCL  05  x  (RCL  03  -  RCL  05)  -  RCL  03 

x  RCL  06  x2)  t  ((RCL  03  -  RCL  05)  x2  +  RCL  06  X2)) 
STO  05  x<>t  STO  06  E 
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APPENDIX  C 


BANDSHAPE  PROGRAM 

i 

This  Fortran  program  plots  Bragg  cell  passband  shapes  for 
arbitrary  integer  ratios  of  interaction  length  (L)  to  charac¬ 
teristic  interaction  length  (Lq) .  Selectable  parameters  are: 

NPTS  =  number  of  plot  points  desired 
*DLTAX  =  normalized  frequency  increment 

/ 

*X(1)  =  normalized  start  frequency  for  plot 

The  number  of  plot  lines  is  controlled  by  the  limit  of  the 
second  DO  loop.  Prospective  users  at  the  Naval  Postgraduate 
School  should  refer  to  the  W.  R.  Chruch  Computer  Center's 
Versatec  plotter  routine  user  manual  for  required  job  control 
cards  and  plot  dimensions. 

*  normalized  to  1 
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C  BAND  SHAPE  —  SINX  SQR  OVER  X  SQR 
DIMENSION  X (0201) ,  P(0201) 
INTEGER* 4  ITB (12) /12*0/ 

REAL *4  RTB (28)/28*0.0/ 
EQUIVALENCE  (TITLE , RTB ( 5 ) ) 

REAL*8  TITLE (12) /'BAND  SHAPE'/ 

NPTX  =201 
DLTAX  =  .01 
PI  =  3.141592654 
X(l)  =  .1 
QPRME  =  2.0  *  PI 
Q  =  0.0 

DO  10  I  =  2,  201 
X(I)  =X(I  -1)  +  DLTAX 
10  CONTINUE 

DO  200  J  =  1,  6 
Q  =  QPRME  +  Q 
DO  100  1=1,  201 
B  =  (Q  *  (X (I)  -  X(I) **2)  x  4.0) 
C  =  SIN (  B  ) 

F (I)  =  (  C/B  ) **2 
100  CONTINUE 

IF  (J  .EQ.  1)  ITB  ( 1 )  =  1 

IF  (J  .EQ.  2)  ITB ( 1)  =  2 

IF  (J  .EQ.  6)  ITB ( 1)  =  3 

ITB (3)  =  8 

ITB  (4)  =  4 

ITB (8)  =  2 

ITB (10)  =  2 

RTB ( 1 )  =  .4 

RTB  ( 2 )  =  .4 

CALL  DRAWP  (NPTS , X ,F , ITB , RTB ) 

200  CONTINUE 
STOP 
END 
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APPENDIX  D 


NARROW  OBSTACLE  DIFFRACTION  PROGRAM 

This  Fortran  program  plots  diffraction  patterns  caused 
by  a  narrow  obstacle  in  a  plane-wave  light  field  using  a 
Simpson's  rule  approximation  algorithm.  Selectable  parameters 
are: 

LMBDA  =  laser  wavelength 

YMIN  =  left  plotting  boundary  limit 

DLTAY  =  obstacle  width 

RZERO  =  distance  from  obstacle  to  image  plane 
TNPP  =  total  number  of  plotting  points 
TNIP  =  total  number  of  integration  points 

Prospective  users  at  the  Naval  Postgraduate  School  should 
refer  to  the  W.  R.  Church  Computer  Center's  Versatec  plotter 
routine  user  manual  for  required  job  control  cards  and  plot 
dimensions . 
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C  NARROW  OBSTACLE  DIFFRACTION  PROGRAM 
DIMENSION  Y(500),  X{500),  X(500) 

INTEGER*4  ITB ( 12 ) /12*0/ 

REAL *4  RTB (28)/28*0.0/ 

EQUIVALENCE  (TITLE , RTB ( 5 ) ) 

REAL*8  TITLE ( 12 )/'DLTAY  EQ  2  MICRONS  RZERO  EQ  19  MM'/ 

REAL  LMBDA 

LMBDA  =  .6328  E-06 

YMIN  =  -40.0  E-06 

DLTAY  =  2.0  E-05 

RZERO  =19.0  E-03 

HPI  =  1.570796327 

TNPP  =  500.0 

TNIP  =  8000.0 

STR  =  SQRT (2.0  /  (LMBDA*RZERO) 

UMIN  =  YMIN  *  STR 
DLTAU  =  DLTAY  *  STR 
DII  =  DLTAU  /  TNIP 
DPP  =  -2.0  *  UMIN  /  TNPP 
1  =  1 

TUMIN  =  UMIN 
PUMIN  =  UMIN 
200  SUM1  =0.0 
SUM2  =0.0 

100  Cl  =  COS (HPI* (TUMIN**2) ) 

C2  =  COS (HPI* (TUMIN+DII) **2) 

C3  =  COS (HPI* (TUMIN+2. 0*DII) **2) 

SUM1  =  (Cl  +  4 . 0*C2  +  C3)*DII/3  +  SUMl 
TUMIN  =  TIMIN  +  2 , 0  *DII 

IF  (TUMIN  .LE.  (PUMIN  +  DLTAU  -  2.0*DII))  GO  TO  100 
TUMIN  =  PUMIN 

300  SI  =  SIN (HPI* (TUMIN* *2) ) 

52  =  SIN (HPI* (TUMIN+DII ) **2) 

53  =  SIN(HPI* (TUMIN+2. 0*DII) **2) 

SUM2  =  (SI  +  4 . 0*S2  +S3 ) *DII/3  +  SUM2 
TUMIN  =  TUMIN  +  2.0*DII 

IF  (TUMIN  .LE.  (PUMIN  +  DLTAU  -  2.0*DII))  GO  TO  300 
IF  (I  .EQ.  251)  GO  TO  999 

X(I)  =  (SUM2-SUM1) **2  +  (2.0  -  SUMl  -SUM2) **2 
Z(I)  =  X(I)  /  4.0 
Z ( 501-1 )  =  Z(I) 

Y ( I)  =  (PUMIN  +  DLTAY  /  2.0)  /  STR 
Y  (501-1)  =  -Y ( I ) 

1  =  1  +  1 

PUMIN  =  PUMIN  +  DPP 

NPTS  =  500 

ITB (3)  =  8 

ITB (4)  =  4 

ITB (8)  =  2 

ITB (10)  =  2 

ITB (11)  =  4 
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RTB(l)  »  1.0  E-05 
RTB ( 2 )  *  .5 

GO  TO  200 

999  CALL  DRAWP  (NPTS,Y,Z,ITB,RTB) 
STOP 
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